Measuring flow rates of multiphase fluids

ABSTRACT

A multiphase fluid is flowed from a flow pipe to a U-bend. Several differential pressures of the multiphase fluid flowing through the flow pipe and U-bend are measured. A total flow rate of the multiphase fluid is determined at least based on the measured differential pressures. In some cases, flow rates of each of the phases of the multiphase fluid can be determined at least based on the measured differential pressures.

TECHNICAL FIELD

This disclosure relates to characterization of the flow of multiphase fluids.

BACKGROUND

A multiphase fluid is a mixture of multiple phases of matter. Multiphase fluids can be non-homogenous and can thus exhibit complex flow characteristics. The characteristics of the flow of multiphase fluids in a conduit, for example, can depend on various factors, such as operating conditions (pressure and temperature), composition (and in turn, physical properties such as density and viscosity) of each of the phases, flow rate of each of the phases, and physical characteristics of the conduit (such as diameter and orientation) through which the multiphase fluid is flowing. Characterization of the flow of multiphase fluids can sometimes be difficult. In some cases, for example, in oil and gas operations, accurate flow metering of multiphase fluid mixtures (such as crude oil, natural gas, and brine) can be important.

SUMMARY

This disclosure describes technologies relating to measuring flow rates of multiphase fluids. Certain aspects of the subject matter described can be implemented as a system. The system includes a U-bend. The U-bend is configured to flow a multiphase fluid. The U-bend includes a first conduit, a second conduit, and a connecting conduit. A cross-sectional flow area of the first conduit is substantially the same as a cross-sectional flow area of the second conduit. The connecting conduit connects the first conduit to the second conduit. The system includes a third conduit that is configured to flow the multiphase fluid. The third conduit has a cross-sectional flow area that is different from the cross-sectional flow area of the first conduit. The system includes a first bend that connects the third conduit to the first conduit. The system includes a fourth conduit that is configured to flow the multiphase fluid. The fourth conduit is parallel to and in-line with the third conduit. The fourth conduit has a cross-sectional flow area that is substantially the same as the cross-sectional flow area of the third conduit. The system includes a second bend that connects the fourth conduit to the second conduit. The system includes a first differential pressure sensor and a second differential pressure sensor. The first differential pressure sensor is coupled to the third conduit and to the first conduit. The first differential pressure sensor is configured to measure a first differential pressure of the multiphase fluid across the first bend. The second differential pressure sensor is coupled to the second conduit and to the fourth conduit. The second differential pressure sensor is configured to measure a second differential pressure of the multiphase fluid across the second bend. The system includes a computer. The computer includes a processor and a computer-readable storage medium. The processor is communicatively coupled to the first differential pressure sensor and to the second differential pressure sensor. The computer-readable storage medium is coupled to the processor and stores programming instructions for execution by the processor. The programming instructions instruct the processor to perform operations. The operations include receiving a first differential pressure signal from the first differential pressure sensor. The first differential pressure signal represents the first differential pressure of the multiphase fluid. The operations include receiving a second differential pressure signal from the second differential pressure sensor. The second differential pressure signal represents the second differential pressure of the multiphase fluid. The operations include receiving, as input, a mixture density of the multiphase fluid. The operations include determining a total flow rate of the multiphase fluid at least based on the first differential pressure, the second differential pressure, and the mixture density of the multiphase fluid.

This, and other aspects, can include one or more of the following features.

In some implementations, the first differential pressure sensor is coupled to the third conduit at a first location and to the first conduit at a second location. In some implementations, the first location and the second location are separated by a specified height with respect to gravity. In some implementations, the second differential pressure sensor is coupled to the second conduit at a third location and to the fourth conduit at a fourth location. In some implementations, the third location and the second location are separated by the specified height with respect to gravity.

In some implementations, the total flow rate of the multiphase fluid is determined as a first total mass flow rate of the multiphase fluid at least based on a first differential pressure equation:

m _(T1) =C _(d1)×ε₁ ×K ₁×√{square root over (2×ρ×(ΔP ₁ −ρ×g×h))},

where m_(T1) is the first total mass flow rate of the multiphase fluid, C_(d1) is a first discharge coefficient, ε₁ is a first expansion factor, K₁ is a first fixed geometry factor proportional to the cross-sectional flow area of the first conduit, ΔP₁ is the first differential pressure, ρ is the mixture density of the multiphase fluid, g is an acceleration due to gravity, and h is the specified height.

In some implementations, he total flow rate of the multiphase fluid is determined as a second total mass flow rate of the multiphase fluid at least based on a second differential pressure equation:

m _(T2) =C _(d2)×ε₂ ×K ₂×√{square root over (2×ρ×(ΔP ₂ −ρ×g×h))},

where m_(T2) is the second total mass flow rate of the multiphase fluid, C_(d2) is a second discharge coefficient, ε₂ is a second expansion factor, K₂ is a second fixed geometry factor proportional to the cross-sectional flow area of the fourth conduit, ΔP₂ is the second differential pressure, ρ is the mixture density of the multiphase fluid, g is the acceleration due to gravity, and h is the specified height.

In some implementations, the operations performed by the processor include comparing the first total mass flow rate of the multiphase fluid and the second total mass flow rate of the multiphase fluid to determine an adjustment of the first discharge coefficient (C_(d1)) or an adjustment of the second discharge coefficient (C_(d2)).

In some implementations, the first total mass flow rate (m_(T1)) is equal to the second total mass flow rate (m_(T2)). In some implementations, the mixture density of the multiphase fluid is a first mixture density of the multiphase fluid. In some implementations, the operations performed by the processor include determining a second mixture density of the multiphase fluid at least based on a combined differential pressure equation:

${\rho_{2} = \frac{{Z \times \Delta P_{1}} - {\Delta P_{2}}}{g \times \left( {1 + Z} \right)}},$

wherein ρ₂ is the second mixture density of the multiphase fluid, ΔP₁ is the first differential pressure, ΔP₂ is the second differential pressure, g is the acceleration due to gravity, and Z is defined by:

${Z = \left( \frac{C_{d1} \times \varepsilon_{1} \times K_{1}}{C_{d2} \times \varepsilon_{2} \times K_{2}} \right)^{2}},$

wherein C_(d1) is the first discharge coefficient, ε₁ is the first expansion factor, K₁ is the first fixed geometry factor, C_(d2) is the second discharge coefficient, ε₂ is the second expansion factor, and K₂ is the second fixed geometry factor.

In some implementations, the first differential pressure sensor is a first circumferential pressure sensor that spans an entire circumference of the first conduit. In some implementations, the second differential pressure sensor is a second circumferential pressure sensor that spans an entire circumference of the second conduit. In some implementations, each of the first bend and the second bend are 90-degree angled bends.

Certain aspects of the subject matter described can be implemented as a method. A multiphase fluid is flowed from a flow pipe to a U-bend. The multiphase fluid includes a gas phase, an aqueous phase, and an oil phase. The U-bend includes a first conduit, a second conduit, and a connecting conduit. The connecting conduit connects the first conduit to the second conduit. The multiphase fluid flowing through the U-bend flows into the first conduit, through the connecting conduit, and out of the second conduit. The multiphase fluid is flowed from the U-bend to the flow pipe. The flow pipe includes a third conduit and a fourth conduit. The third conduit of the flow pipe is connected to the first conduit of the U-bend by a first bend. The fourth conduit of the flow pipe is connected to the second conduit of the U-bend by a second bend. Flowing the multiphase fluid from the flow pipe to the U-bend includes flowing the multiphase fluid from the third conduit (of the flow pipe) through the first bend to the first conduit (of the U-bend). Flowing the multiphase fluid from the U-bend to the flow pipe includes flowing the multiphase fluid from the second conduit (of the U-bend) through the second bend to the fourth conduit (of the flow pipe). A first differential pressure of the multiphase fluid is measured between a first location on the third conduit of the flow pipe and a second location on the first conduit of the U-bend. The first location and the second location are separated by a specified height along the first conduit of the U-bend. A second differential pressure of the multiphase fluid is measured between a third location on the second conduit of the U-bend and a fourth location on the fourth conduit of the flow pipe. The third location and the fourth location are separated by the specified height along the second conduit of the U-bend. A total flow rate of the multiphase fluid is determined at least based on the first differential pressure, the second differential pressure, and a mixture density of the multiphase fluid.

This, and other aspects, can include one or more of the following features.

In some implementations, the total flow rate of the multiphase fluid is determined as a first total mass flow rate of the multiphase fluid at least based on a first differential pressure equation:

m _(T1) =C _(d1)×ε₁ ×K ₁×√{square root over (2×ρ×(ΔP ₁ −ρ×g×h))},

where m_(T1) is the first total mass flow rate of the multiphase fluid, C_(d1) is a first discharge coefficient, ε₁ is a first expansion factor, K₁ is a first fixed geometry factor proportional to the cross-sectional flow area of the first conduit, ΔP₁ is the first differential pressure, ρ is the mixture density of the multiphase fluid, g is an acceleration due to gravity, and h is the specified height.

In some implementations, the total flow rate of the multiphase fluid is determined as a second total mass flow rate of the multiphase fluid at least based on a second differential pressure equation:

m _(T2) =C _(d2)×ε₂ ×K ₂×√{square root over (2×ρ×(ΔP ₂ −ρ×g×h))},

wherein m_(T2) is the second total mass flow rate of the multiphase fluid, C_(d2) is a second discharge coefficient, ε₂ is a second expansion factor, K₂ is a second fixed geometry factor proportional to the cross-sectional flow area of the fourth conduit, ΔP₂ is the second differential pressure, ρ is the mixture density of the multiphase fluid, g is the acceleration due to gravity, and h is the specified height.

In some implementations, the first total mass flow rate (m_(T1)) is equal to the second total mass flow rate (m_(T2)). In some implementations, the mixture density of the multiphase fluid is a first mixture density of the multiphase fluid. In some implementations, the operations performed by the processor include determining a second mixture density of the multiphase fluid at least based on a combined differential pressure equation:

${\rho_{2} = \frac{{Z \times \Delta P_{1}} - {\Delta P_{2}}}{g \times \left( {1 + Z} \right)}},$

where ρ₂ is the second mixture density of the multiphase fluid, ΔP₁ is the first differential pressure, ΔP₂ is the second differential pressure, g is the acceleration due to gravity, and Z is defined by:

${Z = \left( \frac{C_{d1} \times \varepsilon_{1} \times K_{1}}{C_{d2} \times \varepsilon_{2} \times K_{2}} \right)^{2}},$

where C_(d1) is the first discharge coefficient, ε₁ is the first expansion factor, K₁ is the first fixed geometry factor, C_(d2) is the second discharge coefficient, ε₂ is the second expansion factor, and K₂ is the second fixed geometry factor.

In some implementations, the first total mass flow rate of the multiphase fluid and the second total mass flow rate of the multiphase fluid are compared to determine an adjustment of the first discharge coefficient (C_(d1)) or an adjustment of the second discharge coefficient (C_(d2)).

In some implementations, the first differential pressure is measured by a first differential pressure sensor coupled to the third conduit of the flow pipe at the first location and coupled to the first conduit of the U-bend at the second location. In some implementations, the second differential pressure is measured by a second differential pressure sensor coupled to the second conduit of the U-bend at the third location and coupled to the fourth conduit of the flow pipe at the fourth location. In some implementations, the first differential pressure sensor is a first circumferential pressure sensor that spans an entire circumference of the first conduit. In some implementations, the second differential pressure sensor is a second circumferential pressure sensor that spans an entire circumference of the second conduit. In some implementations, each of the first bend and the second bend are 90-degree angled bends.

The details of one or more implementations of the subject matter of this disclosure are set forth in the accompanying drawings and the description. Other features, aspects, and advantages of the subject matter will become apparent from the description, the drawings, and the claims.

DESCRIPTION OF DRAWINGS

FIG. 1A is a schematic diagram of an example apparatus which can be used to measure a density of a multiphase fluid.

FIG. 1B is a schematic diagram of an example apparatus which can be used to measure a density of a multiphase fluid.

FIG. 1C is a schematic diagram of an example apparatus which can be used to measure a density of a multiphase fluid.

FIG. 1D is a schematic diagram of an example apparatus which can be used to measure a density of a multiphase fluid.

FIG. 1E is a schematic diagram of an example apparatus which can be used to measure a flow rate of a multiphase fluid.

FIG. 1F is a schematic diagram of an example apparatus which can be used to measure a flow rate of a multiphase fluid.

FIG. 1G is a schematic diagram of an example apparatus which can be used to measure a water-liquid ratio of a multiphase fluid.

FIG. 1H is a schematic diagram of an example apparatus which can be used to measure a density and a flow rate of a multiphase fluid.

FIG. 1I is a schematic diagram of an example apparatus which can be used to measure a density, a flow rate, and a water-liquid ratio of a multiphase fluid.

FIG. 1J is a schematic diagram of an example apparatus which can be used to measure a density, a flow rate, and a water-liquid ratio of a multiphase fluid.

FIG. 2A is a schematic diagram depicting various example pressure tap arrangements.

FIG. 2B is a schematic diagram depicting various example piping arrangements.

FIG. 2C is a schematic diagram depicting various example piping arrangements.

FIG. 2D is a schematic diagram depicting various example piping arrangements.

FIG. 2E is a schematic diagram depicting various example piping arrangements.

FIG. 2F is a schematic diagram depicting various example sensor arrangements.

FIG. 3A is a flowchart depicting an example multiphase fluid flow model.

FIG. 3B is a flowchart depicting an example multiphase fluid flow model for determining density.

FIG. 3C is a flowchart depicting an example multiphase fluid flow model for determining density.

FIG. 3D is a flowchart depicting an example multiphase fluid flow model for determining density.

FIG. 3E is a flowchart depicting an example multiphase fluid flow model for determining flow rate.

FIG. 3F is a flowchart depicting an example multiphase fluid flow model.

FIG. 3G is a flowchart depicting an example multiphase fluid flow model.

FIG. 4 is a flowchart of an example method for measuring a density of a multiphase fluid.

FIG. 5 is a flowchart of an example method for measuring a flow rate of a multiphase fluid.

FIG. 6 is a flowchart of an example method for measuring a density and a flow rate of a multiphase fluid.

FIG. 7 is a block diagram of an example computer system.

DETAILED DESCRIPTION

This disclosure describes flow characterization of multiphase fluids. The apparatuses described here can be used to accurately measure the total flow rate and the phase flow rates of a multiphase fluid. This information can then be used to determine other properties of the multiphase fluid, for example, based on additional known information. The apparatuses include multiple connected conduits including an elbow (or L-bend) and a U-bend. The apparatuses include multiple pressure sensors for measuring differential pressures across various location pairs of the flowing multiphase fluid. The measured differential pressures can be compared and manipulated to determine the total flow rate and the phase flow rates of the flowing multiphase fluid.

The subject matter described in this disclosure can be implemented in particular implementations, so as to realize one or more of the following advantages. The apparatuses, systems, and methods described here can be implemented independent of the use of radioactive energy sources (such as gamma-ray attenuation), which can be hazardous. The apparatuses, systems, and methods described here can be implemented to accurately determine the total flow rate and the phase flow rates of a multiphase fluid, even in situations in which a gas phase is present in a non-negligible amount. The apparatuses, systems, and methods described here can be implemented to accurately determine the total flow rate and the phase flow rates of a multiphase fluid, while taking into account frictional and shear pressure losses of the multiphase fluid. The apparatuses, systems, and methods described here can be implemented to accurately determine the total flow rate and the phase flow rates of a multiphase fluid, even in complex flow regimes such as slug-flow, plug-flow, and annular flow regimes. The apparatuses, systems, and methods described here can be implemented to accurately determine the total flow rate, the mixture density, and the phase flow rates of a multiphase fluid, even in situations in which there is slip (that is, differential flow velocity) between gas and liquid phases. The apparatuses, systems, and methods described here can take advantage of pipe geometry (and associated fluid dynamics) to measure a relatively simple (and more direct) flow/fluid property (such as differential pressure), rather than a relatively complex (and more indirect) fluid/flow property (such as electromagnetic (spectral) absorption/transmission) to accurately determine the total flow rate and the phase flow rates of a multiphase fluid. The apparatuses, systems, and methods described here can be implemented in conjunction with other flow metering components to supplement and/or enhance the accuracy of the characterization of multiphase fluid flow.

FIG. 1A depicts an apparatus 100A which can be used to measure a density of a multiphase fluid 101. The multiphase fluid 101 can be a 2-phase or a 3-phase fluid. For example, the multiphase fluid 101 can include a gas phase 101 a and a liquid phase 101 b. In some implementations, the liquid phase 101 b includes an aqueous phase 101 b′ and an oil phase 101 b″. For example, the multiphase fluid 101 can include a gas phase 101 a, an aqueous phase 101 b′, and an oil phase 101 b″. As one example, the multiphase fluid 101 can include hydrocarbon gas, hydrocarbon liquid, and water (or brine).

The apparatus 100A includes a U-bend 102 which is configured to flow the multiphase fluid 101. The U-bend 102 includes a first conduit 102 a, a second conduit 102 b, and a connecting conduit 102 c. In some implementations, the second conduit 102 b is parallel to the first conduit 102 a. In some implementations, a longitudinal length of the first conduit 102 a is substantially the same as a longitudinal length of the second conduit 102 b, and a cross-sectional area of the first conduit 102 a is substantially the same as a cross-sectional area of the second conduit 102 b, such that a frictional component of the pressure drop experienced by the multiphase fluid 101 flowing through these components (102 a, 102 b) of the U-bend 102 is substantially the same for each component. For example, each of the first conduit 102 a and the second conduit 102 b have substantially the same inner diameter, d. The connecting conduit 102 c connects the first conduit 102 a to the second conduit 102 b. In some implementations, the first conduit 102 a, the second conduit 102 b, and the connecting conduit 102 c are integrated, such that the U-bend 102 is a singular, unitary body, as opposed to parts that are disjointed and coupled together to form the U-bend 102. In some implementations, the connecting conduit 102 c is perpendicular to the first conduit 102 a and to the second conduit 102 b. The multiphase fluid 101 flowing through the U-bend 102 flows into the first conduit 102 a, through the connecting conduit 102 c, and out of the second conduit 102 b. In some implementations, each of the first conduit 102 a, the second conduit 102 b, and the connecting conduit 102 c are made of the same material, such that the friction experienced by the multiphase fluid 101 flowing through each of the first conduit 102 a, the second conduit 102 b, and the connecting conduit 102 c is substantially the same. In some implementations, the connecting conduit 102 b is connected to the first conduit 102 a by a first curved bend 103 a. In some implementations, the connecting conduit 102 c is connected to the second conduit 102 b by a second curved bend 103 b. In some implementations, the first conduit 102 a, the second conduit 102 b, the connecting conduit 102 c, the first curved bend 103 a, and the second curved bend 103 b are integrated, such that the U-bend 102 is a singular, unitary body, as opposed to parts that are disjointed and coupled together to form the U-bend 102. The first curved bend 103 a and the second curved bend 103 b can be made of the same material as the first conduit 102 a, the second conduit 102 b, and the connecting conduit 102 c, such that the friction experienced by the multiphase fluid 101 flowing through each of the components of the U-bend 102 is substantially the same.

The apparatus 100A includes a first pair of pressure sensors 104 a, 104 b. The pressure sensors 104 a, 104 b are coupled to the first conduit 102 a via pressure ports 105 a, 105 b, respectively. The pressure sensors 104 a, 104 b are cooperatively configured to measure a first differential pressure of the multiphase fluid 101 flowing through the U-bend 102. The pressure ports 105 a, 105 b to which the pressure sensors 104 a, 104 b are coupled, respectively, are separated by a specified differential pressure height, h, along the first tubular 102 a. The apparatus 100A includes a second pair of pressure sensors 104 c, 104 d. The pressure sensors 104 c, 104 d are coupled to the second conduit 102 b via pressure ports 105 c, 105 d. The pressure sensors 104 c, 104 d are cooperatively configured to measure a second differential pressure of the multiphase fluid 101 flowing through the U-bend 102. The pressure ports 105 c, 105 d to which the pressure sensors 104 c, 104 d are coupled, respectively, are separated by the same specified differential pressure height, h, along the second tubular 102 b. In other words, the height difference between the pressure ports 105 a, 105 b of the pressure sensors 104 a, 104 b and the height difference between the pressure ports 105 c, 105 d of the pressure sensors 104 c, 104 d are the same, such that the frictional components of the pressure drops experienced by the multiphase fluid 101 flowing through the first conduit 102 a (measured by pressure sensors 104 a, 104 b) and through the second conduit 102 b (measured by the pressure sensors 104 c, 104 d) are substantially the same. In some implementations, the first pair of pressure sensors 104 a, 104 b is replaced by a first differential pressure sensor that is coupled to the first conduit 102 a via the pressure ports 105 a, 105 b. In some implementations, the second pair of pressure sensors 104 c, 104 d is replaced by a second differential pressure sensor that is coupled to the second conduit 102 b via the pressure ports 105 c, 105 d.

As shown in FIG. 1A, the U-bend 102 can be connected to a flow pipe 150 that flows the multiphase fluid 101. The flow pipe 150 includes a first pipe 150 a and a second pipe 150 b. The terms “pipe” and “conduit” are synonymous. The term “pipe” is being used in relation to the flow pipe 150 and the term “conduit” is being used in relation to the U-bend 102 simply for consistency/clarity and to avoid confusion. In some implementations, the U-bend 102 is connected to the flow pipe 150 by angled bends 153 a, 153 b. For example, the first conduit 102 a of the U-bend 102 is connected to the first pipe 150 a of the flow pipe 150 by the first angled bend 153 a. For example, the second conduit 102 b of the U-bend 102 is connected to the second pipe 150 b of the flow pipe 150 by the second angled bend 153 b. In some implementations, as shown in FIG. 1A, the angled bends 153 a, 153 b are 90-degree angled bends. Although shown as 90-degree angled bends in FIG. 1A, the angled bends 153 a, 153 b can have a different degree angle, such as in a range of from 30 degrees to 150 degrees. In some implementations, the first conduit 102 a, the second conduit 102 b, the connecting conduit 102 c, the first curved bend 103 a, the second curved bend 103 b, the first angled bend 153 a, the second angled bend 153 b, and the flow pipe 150 are integrated, such that the U-bend 102 and the flow pipe 150 form a singular, unitary body, as opposed to parts that are disjointed and coupled together. In some implementations, the angled bends 153 a, 153 b are replaced by curved bends (similar to the curved bends 103 a, 103 b).

In some implementations, the inner diameters (d) of the first conduit 102 a and the second conduit 102 b are different from the inner diameter (D) of the flow pipe 150. A change in cross-sectional flow area between the flow pipe 150 and the U-bend 102 can ensure that the differential pressures of the multiphase fluid 101 flowing across the angled bends 153 a, 153 b are predominantly functions of the total flow rate of the multiphase fluid 101 as opposed to other factors, such as frictional pressure losses. In some implementations, the inner diameters (d) of the first conduit 102 a and the second conduit 102 b can be in a range of from half of the inner diameter (D) of the flow pipe 150 to double the inner diameter (2D) of the flow pipe 150 (that is, D/2≤d≤2D).

In some implementations, each of the pressure sensors (104 a, 104 b, 104 c, 104 d) are at least a threshold distance away from the closest bend. The threshold distance is long enough to allow the flow of the multiphase fluid 101 to fully develop after the disturbance in flow caused by an upstream bend or long enough to allow the flow of the multiphase fluid 101 to remain fully developed before approaching a downstream bend. Without being bound to theory, the method of sensing (for example, the application and/or objective of the measurement) can determine the location at which the sensors should be placed. For the apparatuses and methods described here, the threshold distance can be, for example, four times the inner diameter, d (that is, 4d). For example, if the inner diameter (d) of the U-bend 102 is 1 inch, the threshold distance is 4 inches. For example, the pressure port 105 a, to which the pressure sensor 104 a is coupled, is located on the first conduit 102 a at least 4d distance away from the angled bend 153 a. For example, the pressure port 105 b, to which the pressure sensor 104 b is coupled, is located on the first conduit 102 a at least 4d distance away from the curved bend 103 a. For example, the pressure port 105 c, to which the pressure sensor 104 c is coupled, is located on the second conduit 102 b at least 4d distance away from the curved bend 103 b. For example, the pressure port 105 d, to which the pressure sensor 104 d is coupled, is located on the second conduit 102 b at least 4d distance away from the angled bend 153 b.

As shown in FIG. 1A, the apparatus 100A can be communicatively coupled to a computer 700. The computer 700 includes a processor 705 and a memory 707. The processor 705 can be communicatively coupled to the first pair of pressure sensors 104 a, 104 b (or first differential pressure sensor) and to the second pair of pressure sensors 104 c, 104 d (or second differential pressure sensor). The memory 707 is coupled to the processor 705. The computer 700 is also shown in FIG. 7 and is described in more detail later. The memory 707 stores programming instructions for execution by the processor 705. The programming instructions can instruct the processor 705 to perform various operations. The operations performed by the processor 705 can include receiving a first differential pressure signal from the first pair of pressure sensors 104 a, 104 b (or first differential pressure sensor). The first differential pressure signal can represent the first differential pressure of the multiphase fluid 101 measured by the first pair of pressure sensors 104 a, 104 b (or first differential pressure sensor). The operations performed by the processor 705 can include receiving a second differential pressure signal from the second pair of pressure sensors 104 c, 104 d (or second differential pressure sensor). The second differential pressure signal can represent the second differential pressure of the multiphase fluid 101 measured by the pressure sensors 104 c, 104 d (or second differential pressure sensor). The operations performed by the processor 705 can include determining an overall density of the multiphase fluid 101 at least based on a difference between the first differential pressure and the second differential pressure of the multiphase fluid 101 and the specified differential pressure height, h.

In cases in which the first conduit 102 a and the second conduit 102 b are in a vertical orientation (relative to a reference horizontal, such as ground level), the differential pressure loss of multiphase fluid 101 flowing across a vertical distance can be estimated by Equation 1:

ΔP=ρ×g×h _(v) ×P _(f)  (1)

where ΔP is the total differential pressure loss of the multiphase fluid 101, ρ is the apparent mixture density of the multiphase fluid 101, g is acceleration due to gravity, h_(v), is the difference in vertical height of two reference points (such as the locations of the pressure ports 105 a, 105 b or the locations of the pressure ports 105 c, 105 d), and ΔP_(f) is the pressure loss due to miscellaneous factors, such as frictional losses and multiphase shear. Thus, the first differential pressure measured by the first pair of pressure sensors 104 a, 104 b can be estimated by Equation 1a, and the second differential pressure measured by the second pair of pressure sensors 104 c, 104 d can be estimated by Equation 1b:

ΔP ₁ =ρ×g×h ₁ +ΔP _(f1)  (1a)

ΔP ₂ =−ρ×g×h ₂ +ΔP _(f2)  (1b)

The apparent mixture density of the multiphase fluid (ρ) can be determined by obtaining the difference between the first and second differential pressures (ΔP₁−ΔP₂) and re-arranging the resulting equation to solve for ρ. The difference between the first and second differential pressures is shown by Equation 2a:

ΔP ₁ −ΔP ₂ =ρ×g×h ₁ +ΔP _(f1) +ρ×g×h ₂ −ΔP _(f2)  (2a)

Because the flow characteristics of the first conduit 102 a and the second conduit 102 b are substantially the same, ΔP_(f1) and ΔP_(f2) can be assumed to be substantially the same (ΔP_(f1)≅ΔP_(f2)), thus Equation 2a can be re-written as Equation 2b:

ΔP ₁ −ΔP ₂ =ρ×g×h ₁ +ρ×g×h ₂  (2b)

In cases in which the first differential pressure (ΔP₅₂ ) and the second differential pressure (ΔP_(↓)) are each measured across the same vertical distance, but at opposite flow directions (as is the case for the first pair of pressure sensors 104 a, 104 b and the second pair of pressure sensors 104 c, 104 d shown in FIG. 1A), the first vertical height (h₁) is equal to the second vertical height (h₂)−that is, h₁=h₂. Thus, Equation 2b can be re-written as Equation 2c:

ΔP ₁ −ΔP ₂ =ρ×g×h ₁  (2c)

Inputting h as the vertical height for h₁ and re-arranging Equation 2c to solve for the overall density of the multiphase fluid (ρ) results in Equation 2d:

$\begin{matrix} {\rho = \frac{{\Delta P_{1}} - {\Delta P_{2}}}{2 \times g \times h}} & \left( {2d} \right) \end{matrix}$

The first differential pressure (ΔP₁) is measured by the first pair of pressure sensors 104 a, 104 b (or first differential pressure sensor). The second differential pressure (ΔP₂) is measured by the second pair of pressure sensors 104 c, 104 d (or second differential pressure sensor). The acceleration of gravity is a known universal constant. Thus, ΔP₁, ΔP₂, g, and h can be input to Equation 2d to estimate the mixture density (ρ) of the multiphase fluid 101.

Additionally, information about the frictional/shear pressure losses can be estimated by adding the first differential pressure (ΔP₁) and the second differential pressure (ΔP₂), as shown in Equation 2e:

ΔP ₁ +ΔP ₂ =ρ×g×h ₁ +ΔP _(f1) −ρ×g×h ₂ +ΔP _(f2)  (2e)

As mentioned previously, in cases where the first vertical height (h₁) is equal to the second vertical height (h₂), Equation 2e can be re-written as Equation 2f:

ΔP ₁ +ΔP ₂ =ΔP _(f1) +ΔP _(f2)  (2f)

As mentioned previously, because the flow characteristics of the first conduit 102 a and the second conduit 102 b are substantially the same, ΔP_(f1) and ΔP_(f2) can be assumed to be substantially the same (ΔP_(f1)≅ΔP_(f2)) and can simply be substituted as ΔP_(f), thus Equation 2f can be re-written as Equation 2g:

ΔP₁ +ΔP ₂=2×ΔP _(f)  (2g)

Re-arranging Equation 2g to solve for frictional/shear pressure losses (ΔP_(f)) provides Equation 2h:

$\begin{matrix} {{\Delta P_{f}} = \frac{{\Delta P_{1}} + {\Delta P_{2}}}{2}} & \left( {2h} \right) \end{matrix}$

The information obtained related to frictional/shear pressure losses (ΔP_(f)) can be useful for determining other correlations and flow characteristics of the multiphase fluid 101. For example, if the apparatus 100A is used in conjunction with a Venturi-type flowmeter, the information obtained related to frictional/shear pressure losses (ΔP_(f)) can be used to determine correlations of variables used in the Venturi-type flowmeter calculations. For a Venturi-type flowmeter, Bernoulli's equation can be converted into the form shown in Equation 3a:

m _(T) =C _(d) ×ε×K×√{square root over (2×ρ×(ΔP−k _(f) ×ΔP _(f)))}  (3a)

where m_(T) is the total mass flow rate of the multiphase fluid 101, C_(d) is the discharge coefficient, ε is the expansion factor, K is a fixed geometry factor (K=A/√{square root over (1−β⁴)}, where A is the cross-sectional flow area of the Venturi-type flowmeter's neck (d) and β is the beta ratio (d/D, where d is the neck diameter of the Venturi-type flowmeter and D is the upstream pipe which converges with the Venturi-type flowmeter)), ΔP is measured pressure drop including frictional/shear pressure losses, k_(f) is a geometry-based correction factor to account for difference(s) in the shape/geometry of the Venturi-type flowmeter (for example, as in apparatus 100A), ΔP_(f) is frictional/shear pressure losses (for example, measured by apparatus 100A and Equation 2h), and ρ is the mixture density of the multiphase fluid 101. The overall density (ρ) of the multiphase fluid 101 and frictional/shear pressure losses (ΔP_(f)) determined by the apparatus 100A can, for example, be used in conjunction with the information obtained by the Venturi-type flowmeter (for example, total mass flow rate of the multiphase fluid 101, m_(T)) and be input to Equation 3a.

The operations performed by the processor 705 can include receiving, as input, a gas phase volume fraction (α_(g)) of the multiphase fluid 101. The operations performed by the processor 705 can include receiving, as input, a gas phase density (ρ_(g)) of the gas phase 101 a the multiphase fluid 101. The operations performed by the processor 705 can include receiving, as input, an aqueous phase density (ρ_(a)) of the aqueous phase 101 b′ of the multiphase fluid 101. The operations performed by the processor 705 can include receiving, as input, an oil phase density (ρ_(o)) of the oil phase 101 b″ of the multiphase fluid 101. The operations performed by the processor 705 can include determining an aqueous phase volume fraction (α_(a)) of the multiphase fluid 101 at least based on the gas phase volume fraction (α_(g)), the gas density (ρ_(g)), the aqueous phase density (ρ_(a)), the oil phase density (ρ_(o)), and the overall density (ρ) of the multiphase fluid 101. For example, the gas phase volume fraction (α_(g)), the gas density (ρ_(g)), the aqueous phase density (ρ_(a)), the oil phase density (ρ_(o)), and the overall density (ρ) of the multiphase fluid 101 can be input to Equation 3b and re-arranged to calculate the aqueous phase volume fraction (α_(a)) of the multiphase fluid 101:

ρ=ρ_(g)×α_(g)+ρ_(a)×α_(a)+ρ_(o)×(1−α_(g)−α_(a))  (3b)

The operations performed by the processor 705 can include receiving, as input, an aqueous phase volume fraction (α_(a)) of the multiphase fluid 101. The operations performed by the processor 705 can include determining a gas phase volume fraction (α_(g)) of the multiphase fluid 101 at least based on the aqueous phase volume fraction (α_(a)), the gas density (ρ_(g)), the aqueous phase density (ρ_(a)), the oil phase density (ρ_(o)), and the overall density (ρ) of the multiphase fluid 101. For example, the aqueous phase volume fraction (α_(a)), the gas density (ρ_(g)), the aqueous phase density (ρ_(a)), the oil phase density (ρ_(o)), and the overall density (ρ) of the multiphase fluid 101 can be input to Equation 3b and re-arranged to calculate the gas phase volume fraction (α_(g)) of the multiphase fluid 101.

In some implementations, the first pair of pressure sensors 104 a, 104 b and the second pair of pressure sensors 104 c, 104 d measure the first differential pressure and the second differential pressure, respectively, of the multiphase fluid 101 flowing through the U-bend 102 multiple times across a time duration. Each of the first differential pressures and each of the second differential pressures measured by the first pair of pressure sensors 104 a, 104 b and the second pair of pressure sensors 104 c, 104 d, respectively, can be correlated to time points at which they were measured. Each of the first differential pressures and each of the second differential pressures measured by the first pair of pressure sensors 104 a, 104 b and the second pair of pressure sensors 104 c, 104 d, respectively, can take account for the time delay of the multiphase fluid 101 taking time to travel from the first conduit 102 a (where the first pair of pressure sensors 104 a, 104 b are located), through the connecting conduit 102 c, and to the second conduit 102 b (where the second pair of pressure sensors 104 c, 104 d are located). For example, the measured differential pressures can be time-corrected to account for the delay in the multiphase fluid 101 traveling through the U-bend 102. The time-series data can, for example, be stored in the memory 707. The time-series data can, for example, be used to train and/or build a neural network-based classification model, which can be used to accurately identify flow regimes (for example, bubble flow, mist flow, slug flow, churn flow, annular flow, stratified flow, or intermittent flow). For example, the time-series data can be converted into a spectrogram (for example, by using a Morlet wavelet transform), and a single-layer two-dimensional image can be compiled from the spectrogram. A neural network (machine learning model) can then analyze the image and identify the flow regime of the multiphase fluid 101 flowing through the U-bend 102, for example, based on comparison to historical data. The time-series data can, for example, be used to train and/or build a neural network-based classification model, which can be used to estimate bulk flow velocity of the multiphase fluid 101 flowing through the U-bend 102. For example, the time-series data can be converted into a spectrogram (for example, by using a Morlet wavelet transform), and a single-layer two-dimensional image can be compiled from the spectrogram. A neural network (machine learning model) can then analyze the image and estimate the bulk flow velocity and/or flow rate of the multiphase fluid 101 flowing through the U-bend 102, for example, based on comparison to historical data.

FIG. 1B depicts an apparatus 100B which can be used to measure a density of a multiphase fluid 101. The apparatus 100B can be substantially similar to the apparatus 100A shown in FIG. 1A. In comparison, the apparatus 100B omits the connecting conduit (102 c shown in FIG. 1A). The apparatus 100B includes the U-bend 102, which includes the first conduit 102 a and the second conduit 102 b. The first conduit 102 a is connected to the second conduit 102 b by a curved bend 103. In some implementations, the first conduit 102 a, the second conduit 102 b, and the curved bend 103 are integrated, such that the U-bend 102 is a singular, unitary body, as opposed to parts that are disjointed and coupled together to form the U-bend 102. The curved bend 103 can be made of the same material as the first conduit 102 a and the second conduit 102 b, such that the friction experienced by the multiphase fluid 101 flowing through each of the components of the U-bend 102 is substantially the same.

The apparatus 100B can include a first differential pressure sensor 108 _(a) coupled to the first conduit 102 a via pressure ports 105 a and 105 b. The first differential pressure sensor 108 _(a)is configured to measure a first differential pressure of the multiphase fluid 101 flowing through the U-bend 102. The pressure ports 105 a and 105 b are separated by the specified differential pressure height, h, with respect to gravity. The apparatus 100B can include a second differential pressure sensor 108 _(b) coupled to the second conduit 102 b via pressure ports 105 c and 105 d. The second differential pressure sensor 108 _(b) is configured to measure a second differential pressure of the multiphase fluid 101 flowing through the U-bend 102. The pressure ports 105 c and 105 d are separated by the specified differential pressure height, h, with respect to gravity. As shown in FIG. 1B, the pressure port 105 a does not necessarily need to be at the same vertical height as pressure port 105 d, and the pressure port 105 b does not necessarily need to be at the same vertical height as pressure port 105 c. Instead, the vertical distance between pressure port 105 a and pressure port 105 b should match the vertical distance between pressure port 105 c and pressure port 105 d. The height difference between the pressure ports 105 a, 105 b of the first differential pressure sensor 108 _(a) and the height difference between the pressure ports 105 c, 105 d of the second differential pressure sensor 108 _(b) are the same, such that the frictional components of the pressure drops experienced by the multiphase fluid 101 flowing through the first conduit 102 a (measured by the first differential pressure sensor 108 _(a)) and through the second conduit 102 b (measured by the second differential pressure sensor 108 _(b)) are substantially the same.

Although not shown in FIG. 1B, the apparatus 100B can be coupled to a flow pipe (such as the flow pipe 150 shown in FIG. 1A), for example, by angled or curved bends. Although not shown in FIG. 1B, the first and second differential pressure sensors 108 _(a), 108 _(b) can be coupled to a computer (such as the computer 700), to determine characteristics (such as a mixture density, a total mass flow rate, an aqueous phase volume fraction, a gas phase volume fraction, or any combination of these) of the multiphase fluid 101 flowing through the U-bend 102 (for example, using Equations 1, 1a, 1b, 2a, 2b, 2c, 2d, 2e, 2f, 2g, 2h, 3a, 3b, or any combination of these).

FIG. 1C depicts an apparatus 100C which can be used to measure a density of a multiphase fluid 101. The apparatus 100C can be substantially similar to the apparatus 100A shown in FIG. 1A. In comparison, the apparatus 100C includes multiple pairs of pressure sensors coupled to the first conduit 102 a and multiple pairs of pressure sensors coupled to the second conduit 102 b. The apparatus 100C shown in FIG. 1C includes two pairs of pressure sensors (104 a, 104 b and 104 a′, 104 b′) coupled to the first conduit 102 a installed in a parallel configuration, but can, in other embodiments, include additional pairs of pressure sensors. The apparatus 100C shown in FIG. 1C includes two pairs of pressure sensors (104 c, 104 d and 104 c′, 104 d′) coupled to the second conduit 102 b installed in a parallel configuration, but can, in other embodiments, include additional pairs of pressure sensors. The redundancy in pressure sensors can be useful in providing more accurate pressure measurements (for example, by averaging the measured pressures). The redundancy in pressure sensors can also be useful in identifying malfunctions in one or more of the pressure sensors, for example, if one or more of the pressure sensors provide pressure readings that deviate from the remaining pressure sensors by more than a threshold deviation. In some implementations, each pair of pressure sensors can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 a, 104 b can be replaced by a first differential pressure sensor, the pair of pressure sensors 104 a′, 104 b′ can be replaced by a second differential pressure sensor, the pair of pressure sensors 104 c, 104 d can be replaced by a third differential pressure sensor, and the pair of pressure sensors 104 c′, 104 d′ can be replaced by a fourth differential pressure sensor. The redundancy in differential pressure sensors can be useful in providing a wider range of measurements (for example, by using sensors with different, but overlapping ranges).

FIG. 1D depicts an apparatus 100D which can be used to measure a density of a multiphase fluid 101. The apparatus 100D can be substantially similar to the apparatus 100A shown in FIG. 1A. In comparison, the apparatus 100D includes more than two pressure sensors coupled to the first conduit 102 a and more than two pressure sensors coupled to the second conduit 102 b. Although shown in FIG. 1D as including six pressure sensors (104 a, 104 b, 104 c, 104 d, 104 e, 104 f) coupled to the first conduit 102 a, the apparatus 100D can include fewer pressure sensors (such as three, four, or five) or more pressure sensors (such as seven or more than seven) coupled to the second conduit 102 b. Although shown in FIG. 1D as including six pressure sensors (104 a′, 104 b′, 104 c′, 104 d′, 104 e′, 104 f′) coupled to the second conduit 102 b, the apparatus 100D can include fewer pressure sensors (such as three, four, or five) or more pressure sensors (such as seven or more than seven) coupled to the second conduit 102 b. Ideally, the apparatus 100D would include the same number of pressure sensors coupled to the first conduit 102 a as coupled to the second conduit 102 b. As shown in FIG. 1D, the pressure sensors (104 a, 104 b, 104 c, 104 d, 104 e, 104 f) are installed on the first conduit 102 a in a series configuration, and the pressure sensors (104 a′, 104 b′, 104 c′, 104 d′, 104 e′, 104 f′) are installed on the second conduit 102 b in a series configuration.

In some cases, the pressure sensors 104 a and 104 b can be considered a first pair for the first conduit 102 a, corresponding to the pressure sensors 104 a′ and 104 b′ as a corresponding first pair for the second conduit 102 b. In some cases, the pressure sensors 104 b and 104 c can be considered a second pair for the first conduit 102 a, corresponding to the pressure sensors 104 b′ and 104 c′ as a corresponding second pair for the second conduit 102 b. In some cases, the pressure sensors 104 c and 104 d can be considered a third pair for the first conduit 102 a, corresponding to the pressure sensors 104 c′ and 104 d′ as a corresponding third pair for the second conduit 102 b. In some cases, the pressure sensors 104 d and 104 e can be considered a fourth pair for the first conduit 102 a, corresponding to the pressure sensors 104 d′ and 104 e′ as a corresponding fourth pair for the second conduit 102 b. In some cases, the pressure sensors 104 e and 104 f can be considered a fifth pair for the first conduit 102 a, corresponding to the pressure sensors 104 e′ and 104 f′ as a corresponding fifth pair for the second conduit 102 b. These pairs can be used to measure a piecewise pressure gradient along the first and second conduits 102 a, 102 b. This piecewise pressure gradient data can, for example, be averaged to more accurately calculate the density of the multiphase fluid 101 flowing through the U-bend.

In some cases, the pressure sensors 104 a and 104 c can be considered a first pair for the first conduit 102 a, corresponding to the pressure sensors 104 a′ and 104 c′ as a corresponding first pair for the second conduit 102 b. In some cases, the pressure sensors 104 b and 104 d can be considered a second pair for the first conduit 102 a, corresponding to the pressure sensors 104 b′ and 104 d′ as a corresponding second pair for the second conduit 102 b. In some cases, the pressure sensors 104 c and 104 e can be considered a third pair for the first conduit 102 a, corresponding to the pressure sensors 104 c′ and 104 e′ as a corresponding third pair for the second conduit 102 b. In some cases, the pressure sensors 104 d and 104 f can be considered a fourth pair for the first conduit 102 a, corresponding to the pressure sensors 104 d′ and 104 f′ as a corresponding fourth pair for the second conduit 102 b. Similarly, these pairs can be used to measure a piecewise pressure gradient along the first and second conduits 102 a, 102 b. This piecewise pressure gradient data can, for example, be averaged to more accurately calculate the density of the multiphase fluid 101 flowing through the U-bend.

In some cases, the pressure sensors 104 a and 104 d can be considered a first pair for the first conduit 102 a, corresponding to the pressure sensors 104 a′ and 104 d′ as a corresponding first pair for the second conduit 102 b. In some cases, the pressure sensors 104 b and 104 e can be considered a second pair for the first conduit 102 a, corresponding to the pressure sensors 104 b′ and 104 e′ as a corresponding second pair for the second conduit 102 b. In some cases, the pressure sensors 104 c and 104 f can be considered a third pair for the first conduit 102 a, corresponding to the pressure sensors 104 c′ and 104 f′ as a corresponding third pair for the second conduit 102 b. Similarly, these pairs can be used to measure a piecewise pressure gradient along the first and second conduits 102 a, 102 b. This piecewise pressure gradient data can, for example, be averaged to more accurately calculate the density of the multiphase fluid 101 flowing through the U-bend.

In some cases, the pressure sensors 104 a and 104 e can be considered a first pair for the first conduit 102 a, corresponding to the pressure sensors 104 a′ and 104 e′ as a corresponding first pair for the second conduit 102 b. In some cases, the pressure sensors 104 b and 104 f can be considered a second pair for the first conduit 102 a, corresponding to the pressure sensors 104 b′ and 104 f′ as a corresponding second pair for the second conduit 102 b. Similarly, these pairs can be used to measure a piecewise pressure gradient along the first and second conduits 102 a, 102 b. This piecewise pressure gradient data can, for example, be averaged to more accurately calculate the density of the multiphase fluid 101 flowing through the U-bend.

In some cases, the pressure sensors 104 a and 104 f can be considered a first pair for the first conduit 102 a, corresponding to the pressure sensors 104 a′ and 104 f′ as a corresponding first pair for the second conduit 102 b. In such cases, the apparatus 100D is substantially the same as the apparatus 100A shown in FIG. 1A. Any combinations or sub-combinations of the aforementioned pairs or pressure sensors can be considered as suitable implementations of the apparatus 100D.

In some implementations, any pair of pressure sensors coupled to the first tubular 102 a can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 a, 104 b can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 a, 104 c can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 a, 104 d can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 a, 104 e can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 a, 104 f can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 b, 104 c can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 b, 104 d can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 b, 104 e can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 b, 104 f can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 c, 104 d can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 c, 104 e can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 c, 104 f can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 d, 104 e can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 d, 104 f can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 e, 104 f can be replaced by a differential pressure sensor.

In some implementations, any pair of pressure sensors coupled to the second tubular 102 b can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 a′, 104 b′ can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 a′, 104 c′ can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 a′, 104 d′ can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 a′, 104 e′ can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 a′, 104 f′ can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 b′, 104 c′ can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 b′, 104 d′ can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 b′, 104 e′ can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 b′, 104 f′ can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 c′, 104 d′ can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 c′, 104 e′ can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 c′, 104 f′ can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 d′, 104 e′ can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 d′, 104 f′ can be replaced by a differential pressure sensor. For example, the pair of pressure sensors 104 e′, 104 f′ can be replaced by a differential pressure sensor.

FIG. 1E depicts an apparatus 100E which can be used to measure a total mass flow rate of a multiphase fluid 101. The apparatus 100E can be substantially similar to the apparatus 100A shown in FIG. 1A. The first pair of pressure sensors 104 a, 104 b and the second pair of pressure sensors 104 c, 104 d are located at different locations on the apparatus 100E in comparison to the apparatus 100A. The first pressure sensor 104 a is coupled to the pressure port 105 a on the first pipe 150 a (third conduit) of the flow pipe 150. The second pressure sensor 104 b is coupled to the pressure port 105 b on the first conduit 102 a of the U-bend 102. The third pressure sensor 104 c is coupled to the pressure port 105 c on the second conduit 102 b of the U-bend 102. The fourth pressure sensor 104 d is coupled to pressure port 105 d on the the second pipe 150 b (fourth conduit) of the flow pipe 150. In some implementations, the first pair of pressure sensors 104 a, 104 b is replaced by a first differential pressure sensor. In some implementations, the second pair of pressure sensors 104 c, 104 d is replaced by a second differential pressure sensor.

The first pair of pressure sensors 104 a, 104 b (or first differential pressure sensor) are configured to measure a first differential pressure of the multiphase fluid 101. The first differential pressure includes the pressure losses of the multiphase fluid 101 flowing through the first angled bend 153 a. The second pair of pressure sensors 104 c, 104 d (or second differential pressure sensor) are configured to measure a second differential pressure of the multiphase fluid 101. The second differential pressure includes the pressure losses of the multiphase fluid 101 flowing through the second angled bend 153 b. The operations performed by the processor 705 can include receiving a first differential pressure signal from the first pair of pressure sensors 104 a, 104 b (or first differential pressure sensor). The first differential pressure signal can represent the first differential pressure of the multiphase fluid 101 measured by the pressure sensors 104 a, 104 b (or first differential pressure sensor). The operations performed by the processor 705 can include receiving a second differential pressure signal from the second pair of pressure sensors 104 c, 104 d (or second differential pressure sensor). The second differential pressure signal can represent the second differential pressure of the multiphase fluid 101 measured by the pressure sensors 104 c, 104 d (or second differential pressure sensor). The operations performed by the processor 705 can include determining an overall flow rate of the multiphase fluid 101 at least based on the first differential pressure, the second differential pressure, a difference between the first differential pressure and the second differential pressure of the multiphase fluid 101 and the overall density of the multiphase fluid 101. The overall density of the multiphase fluid 101 can be measured, for example, by a separate multiphase flowmeter. For example, the apparatus 100E can include features from the apparatus 100A, which can be used to determine the mixture density (ρ) of the multiphase fluid 101.

The overall flow rate of the multiphase fluid 101 can be determined in various ways. As a first example, the overall flow rate of the multiphase fluid 101 (ignoring frictional pressure loss) can be determined using the first differential pressure by Equation 4a:

m _(T1) −C _(d1)×ε₁ ×K ₁×√{square root over (2×ρ×(ΔP ₁ ρ×g×h _(v)))}  (4a)

where m_(T1) is the total mass flow rate of the multiphase fluid 101, C_(d1) is an empirically characterized discharge coefficient, ε₁ is an expansion factor, K₁ is a fixed geometry factor (K₁=A₁/√{square root over (1−β₁ ⁴)}, where A₁ is the cross-sectional flow area of the throat of the first angled bend 153 a and β₁ is the beta ratio (d/D, where d is the inner diameter of the first conduit 102 a and D is the inner diameter of the first pipe 150 a)), ΔP₁ is the first differential pressure measured by the pressure sensors 104 a, 104 b (or the first differential pressure sensor), ρ is the mixture density of the multiphase fluid 101, g is the acceleration due to gravity, and h_(v) is the vertical height difference between the pressure ports 105 a, 105 b to which the pressure sensors 104 a, 104 b (or the first differential pressure sensor) are coupled. The discharge coefficient (C_(d1)) can be determined empirically, for example, based on historical data or correcting with real-time data.

As a second example, the total mass flow rate of the multiphase fluid 101 (ignoring frictional pressure loss) can be determined using the second differential pressure by Equation 4b:

m _(T2) =C _(d2)×ε₂ ×K ₂×√{square root over (2×ρ×(ΔP ₂ ρ×g×h _(v)))}  (4b)

where m_(T2) is the total mass flow rate of the multiphase fluid 101, C_(d2) is an empirically characterized discharge coefficient, ε₂ is an expansion factor, K₂ is a fixed geometry factor (K₂=A₂/√{square root over (1−β₂ ⁴)}, where A₂ is the cross-sectional flow area of the throat of the second angled bend 153 b and β₂ is the beta ratio (D/d, where d is the inner diameter of the second conduit 102 b and D is the inner diameter of the second pipe 150 b)), ΔP₂ is the second differential pressure measured by the pressure sensors 104 c, 104 d (or the second differential pressure sensor), ρ is the mixture density of the multiphase fluid 101, g is the acceleration due to gravity, and h_(v) is the vertical height difference between the pressure ports 105 c, 105 d to which the pressure sensors 104 c, 104 d (or the second differential pressure sensor) are coupled. The discharge coefficient (C_(d2)) can be determined empirically, for example, based on historical data or correcting with real-time data.

The operations performed by the processor 705 can include comparing the total mass flow rates (m_(T1), m_(T2)) determined by Equations 4a and 4b to determine an adjustment of the discharge coefficients (C_(d1), C_(d2)). Equation 4a can be re-written as: m_(T1)=C_(d1)×Y₁, where Y₁ simply substitutes the right-hand side of Equation 4a (including measured parameters ΔP₁ and ρ and predetermined parameters ε₁ and K₁) excluding the discharge coefficient C_(d1). Equation 4b can be re-written as: m_(T2)=C_(d2)×Y₂, where Y₂ simply substitutes the right-hand side of Equation 4b (including measured parameters ΔP₂ and ρ and predetermined parameters ε₂ and K₂) excluding the discharge coefficient C_(d1). Because the total mass flow rate does not change as the multiphase fluid 101 flows through the apparatus 100E, m_(T1)=m_(T2). Thus, C_(d1)×Y₁=C_(d2)×Y₂. Because of this relationship, given C_(d1) and the measured parameters Y₁ and Y₂, C_(d2) can be calculated. Similarly, given C_(d2) and the measured parameters Y₁ and Y₂, C_(d1) can be calculated. As such, the discharge coefficients (C_(d1), C_(d2)) can be verified for consistency and accuracy, as an in situ self-check. Other techniques (such as least-squares minimization) can be additionally or alternatively used to accurately calculate the discharge coefficients (C_(d1), C_(d2)).

The first differential pressure (ΔP₁) measured by the pressure sensors 104 a, 104 b (or the first differential pressure sensor) and the second differential pressure (ΔP₂) measured by the pressure sensors 104 c, 104 d (or the second differential pressure sensor) can be used to calculate the mixture density (ρ) of the multiphase fluid 101. Equation 4a can be re-arranged as Equation 5a:

$\begin{matrix} {{{\Delta P_{1}} - {\rho \times g \times h_{v}}} = \frac{m_{T1}^{2}}{2 \times \rho \times \left( {C_{d1} \times \varepsilon_{1} \times K_{1}} \right)^{2}}} & \left( {5a} \right) \end{matrix}$

Equation 4b can be re-arranged as Equation 5b:

$\begin{matrix} {{{\Delta P_{2}} + {\rho \times g \times h_{v}}} = \frac{m_{T2}^{2}}{2 \times \rho \times \left( {C_{d2} \times \varepsilon_{2} \times K_{2}} \right)^{2}}} & \left( {5b} \right) \end{matrix}$

As mentioned previously, the total mass flow rate does not change as the multiphase fluid 101 flows through the apparatus 100E, so m_(T1)=m_(T2). For simplicity, the variable Z is defined by Equation 5c:

$\begin{matrix} {Z = \left( \frac{C_{d1} \times \varepsilon_{1} \times K_{1}}{C_{d2} \times \varepsilon_{2} \times K_{2}} \right)^{2}} & \left( {5c} \right) \end{matrix}$

Combining Equations 5a, 5b, and 5c and then solving for the mixture density (ρ) results in Equation 5d:

$\begin{matrix} {\rho = \frac{{Z \times \Delta P_{1}} - {\Delta P_{2}}}{g \times \left( {1 + Z} \right)}} & \left( {5d} \right) \end{matrix}$

Equation 5d can be used, for example, to verify the accuracy of the mixture density (ρ) determined by any one of apparatuses 100A, 100B, 100C, or 100D.

In some implementations, a distance between the first pressure sensor 104 a and the second pressure sensor 104 b is substantially the same as a distance between the third pressure sensor 104 c and the fourth pressure sensor 104 d. In a specific example, the apparatus 100E has the following dimensions. The inner diameter (d) of the U-bend 102 (the first conduit 102 a, the second conduit 102 b, and the connecting conduit 102 c) is about 2.067 inches (nominal pipe size (NPS) 2″ Schedule 40). The inner diameter (D) of the first pipe 150 a and the second pipe 150 b is about 3.068″ (NPS 3″ Schedule 40). The dimension A is the length of straight pipe (that is, pipe that is free of any bends) of the first pipe 150 a upstream of the first pressure sensor 104 a. The dimension A is also the length of straight pipe of the second pipe 150 b downstream of the fourth pressure sensor 104 d. The dimension A is a minimum length of straight pipe that ensures minimal interference of turbulence and friction which may affect the measurements taken by the first pair of pressure sensors 104 a, 104 b and the measurements taken by the second pair of pressure sensors 104 c, 104 d. In some implementations, the dimension A is at least five times the inner diameter (D) of the first pipe 150 a (A≥5D). The dimension B is the length of straight pipe of the first pipe 150 a downstream of the first pressure sensor 104 a and upstream of the first angled bend 153 a. The dimension B is also the length of straight pipe of the second pipe 150 b downstream of the second angled bend 153 b and upstream of the fourth pressure sensor 104 d. Similarly, the dimension B is a minimum length of straight pipe that ensures minimal interference of turbulence and friction which may affect the measurements taken by the first pair of pressure sensors 104 a, 104 b and the measurements taken by the second pair of pressure sensors 104 c, 104 d. In some implementations, the dimension B is at least five times the inner diameter (D) of the first pipe 150 a (B≥5D). The dimension C is the length of straight pipe of the first conduit 102 a downstream of the first angled bend 153 a and upstream of the second pressure sensor 104 b. The dimension C is also the length of straight pipe of the second conduit 102 b downstream of the third pressure sensor 104 c and upstream of the second angled bend 153 b. Similarly, the dimension C is a minimum length of straight pipe that ensures minimal interference of turbulence and friction which may affect the measurements taken by the first pair of pressure sensors 104 a, 104 b and the measurements taken by the second pair of pressure sensors 104 c, 104 d. In some implementations, the dimension C is at least five times the inner diameter (d) of the U-bend 102 (C≥5d).

FIG. 1F depicts an apparatus 100F which can be used to measure a total mass flow rate of a multiphase fluid 101. The apparatus 100F can be substantially similar to the apparatus 100E shown in FIG. 1E. In comparison, the apparatus 100F omits the connecting conduit (102 c shown in FIG. 1E). The apparatus 100F includes the U-bend 102, which includes the first conduit 102 a and the second conduit 102 b. The first conduit 102 a is connected to the second conduit 102 bby a curved bend 103. In some implementations, the first conduit 102 a, the second conduit 102 b, and the curved bend 103 are integrated, such that the U-bend 102 is a singular, unitary body, as opposed to parts that are disjointed and coupled together to form the U-bend 102. The curved bend 103 can be made of the same material as the first conduit 102 a and the second conduit 102 b, such that the friction experienced by the multiphase fluid 101 flowing through each of the components of the U-bend 102 is substantially the same.

As shown in FIG. 1F, the U-bend 102 can be connected to a flow pipe 150 that flows the multiphase fluid 101. The flow pipe 150 includes the first pipe 150 a and the second pipe 150 b. In some implementations, the U-bend 102 is connected to the flow pipe 150 by curved bends 153 c, 153 d. For example, the first conduit 102 a of the U-bend 102 is connected to the first pipe 150 a of the flow pipe 150 by the first curved bend 153 c. For example, the second conduit 102 b of the U-bend 102 is connected to the second pipe 150 b of the flow pipe 150 by the second curved bend 153 d. In some implementations, the first conduit 102 a, the second conduit 102 b, the curved bend 103, the first curved bend 153 c, the second curved bend 153 d, and the flow pipe 150 are integrated, such that the U-bend 102 and the flow pipe 150 form a singular, unitary body, as opposed to parts that are disjointed and coupled together.

The apparatus 100F can include a first differential pressure sensor 108 a coupled to the first pipe 150 a and to the first conduit 102 a via pressure ports 105 a and 105 b, respectively. The first differential pressure sensor 108 a is configured to measure a first differential pressure of the multiphase fluid 101 flowing across the first curved bend 153 c. The pressure ports 105 a and 105 b are separated by the specified differential pressure height, h, with respect to gravity. The apparatus 100B can include a second differential pressure sensor 108 b coupled to the second conduit 102 b and to the second pipe 150 b via pressure ports 105 c and 105 d, respectively. The second differential pressure sensor 108 b is configured to measure a second differential pressure of the multiphase fluid 101 flowing across the second curved bend 153 d. The pressure ports 105 c and 105 d are separated by the specified differential pressure height, h, with respect to gravity. As shown in FIG. 1F, the pressure port 105 a does not necessarily need to be at the same vertical height as pressure port 105 d, and the pressure port 105 b does not necessarily need to be at the same vertical height as pressure port 105 c. Instead, the vertical distance between pressure port 105 a and pressure port 105 b should match the vertical distance between pressure port 105 c and pressure port 105 d. The height difference between the pressure ports 105 a, 105 b of the first differential pressure sensor 108 a and the height difference between the pressure ports 105 c, 105 d of the second differential pressure sensor 108 b are the same, and the horizontal distance between the pressure ports 105 a, 105 b of the first differential pressure sensor 108 a and the horizontal distance between the pressure ports 105 c, 105 d of the second differential pressure sensor 108 b are the same, such that the frictional components of the pressure drops experienced by the multiphase fluid 101 flowing across the first curved bend 153 c (measured by the first differential pressure sensor 108 a) and through the second curved bend 153 d (measured by the second differential pressure sensor 108 b) are substantially the same.

Although not shown in FIG. 1F, the first and second differential pressure sensors 108 a, 108 b can be coupled to a computer (such as the computer 700), to determine characteristics (such as a mixture density, a total mass flow rate, or both) of the multiphase fluid 101 flowing through the U-bend 102 (for example, using Equations 4a, 4b, 5a, 5b, 5c, 5d, or any combination of these).

FIG. 1G depicts an apparatus 100G which can be used to measure a water-liquid ratio of a multiphase fluid 101. The apparatus 100G can be substantially similar to the apparatus 100A shown in FIG. 1A. In comparison to the apparatus 100A, the apparatus 100G includes dynamic pressure sensors 106 a, 106 b, 106 c that are coupled to the U-bend 102. Although shown in FIG. 1G as including three dynamic pressure sensors 106 a, 106 b, 106 c, the apparatus 100G can include fewer dynamic pressure sensors (such as two dynamic pressure sensors) or more dynamic pressure sensors (such as four dynamic pressure sensors or more than four dynamic pressure sensors). The dynamic pressure sensors 106 a, 106 b, 106 c measure a dynamic pressure of the multiphase fluid 101 flowing through the U-bend 102, as opposed to a static pressure of the multiphase fluid 101. The dynamic pressure sensors 106 a, 106 b, 106 c can be located at any locations along the U-bend 102. In some implementations, the dynamic pressure sensors 106 a, 106 b, 106 c are coupled to the first conduit 102 a of the U-bend 102. In some implementations, the dynamic pressure sensors 106 a, 106 b, 106 c are coupled to the second conduit 102 b of the U-bend 102. In some implementations, a first portion of the dynamic pressure sensors 106 a, 106 b, 106 c (for example, one or two of the dynamic pressure sensors 106 a, 106 b, 106 c) is coupled to the first conduit 102 a of the U-bend 102, and a remaining portion of the dynamic pressure sensors 106 a, 106 b, 106 c is coupled to the second conduit 102 b of the U-bend 102. The dynamic pressure sensors 106 a, 106 b, 106 c can be evenly or randomly distributed along the U-bend 102.

The dynamic pressure measurements taken by the dynamic pressure sensors 106 a, 106 b, 106 c can be cross-correlated (for example, by the processor 705) to determine a speed of sound of the multiphase fluid 101 flowing through the U-bend 102. The speed of sound of the multiphase fluid 101 and the overall density (ρ) of the multiphase fluid 101 can then be used to determine the water-liquid ratio of the multiphase fluid 101. The water-liquid ratio of the multiphase fluid 101 is a ratio of a volume of the aqueous phase 101 b′ to a sum of the volume of the aqueous phase 101 b′ and a volume of the oil phase 101 b″ (liquid phase 101 b) (that is, the ratio of (volume of 101 b′)/(volume of 101 b′+volume of 101 b″), which is the same as the ratio of (volume of 101 b′)/(volume of 101 b)). For example, the speed of sound of the multiphase fluid 101 and the overall density (ρ) of the multiphase fluid 101 can be compared to a plotted graph of various phase curves (gas/oil, gas/aqueous, oil/aqueous) to determine the water-liquid ratio.

FIG. 1H depicts an apparatus 100H which can be used to measure a density, a total flow rate, and phase flow rates of a multiphase fluid 101. The apparatus 100H can be substantially similar to the apparatus 100A shown in FIG. 1A and the apparatus 100E shown in FIG. 1E. For example, the apparatus 100H can be a combination of the apparatus 100A and the apparatus 100E. The first pressure sensor 104 a is coupled to a first pressure port 105 a on the first pipe 150 a (third conduit) of the flow pipe 150. The second pressure sensor 104 b is coupled to a second pressure port 105 b on the first conduit 102 a of the U-bend 102. The third pressure sensor 104 c is coupled to a third pressure port 105 c on the second conduit 102 b of the U-bend 102. The fourth pressure sensor 104 d is coupled to a fourth pressure port 105 d on the second pipe 150 b (fourth conduit) of the flow pipe 150. The fifth and sixth pressure sensors 104 e, 104 f are coupled to fifth and sixth pressure ports 105 e, 105 f, respectively, on the first conduit 102 a of the U-bend 102. The fifth and sixth pressure ports 105 e, 105 f are separated by the specified height, h, along the first conduit 102 a. The seventh and eighth pressure sensors 104 g, 104 h are coupled to seventh and eighth pressure ports 105 g, 105 h, respectively, on the second conduit 102 b of the U-bend 102. The seventh and eighth pressure ports 105 g, 105 h are separated by the specified height, h, along the second conduit 102 b. In some implementations, the first and second pressure sensors 104 a, 104 b are replaced by a first differential pressure sensor. In some implementations, the third and fourth pressure sensors 104 c, 104 d are replaced by a second differential pressure sensor. In some implementations, the fifth and sixth pressure sensors 104 e, 104 f are replaced by a third differential pressure sensor. In some implementations, the seventh and eighth pressure sensors 104 g, 104 h are replaced by a fourth differential pressure sensor.

The first pair of pressure sensors 104 a, 104 b (or first differential pressure sensor) are configured to measure a first differential pressure of the multiphase fluid 101. The first differential pressure includes the pressure losses of the multiphase fluid 101 flowing through the first angled bend 153 a. The second pair of pressure sensors 104 c, 104 d (or second differential pressure sensor) are configured to measure a second differential pressure of the multiphase fluid 101. The second differential pressure includes the pressure losses of the multiphase fluid 101 flowing through the second angled bend 153 b. The third pair of pressure sensors 104 e, 104 f (or third differential pressure sensor) are configured to measure a third differential pressure of the multiphase fluid 101. The fourth pair of pressure sensors 104 g, 104 h (or fourth differential pressure sensor) are configured to measure a fourth differential pressure of the multiphase fluid 101. The operations performed by the processor 705 can include receiving a first differential pressure signal from the first pair of pressure sensors 104 a, 104 b (or first differential pressure sensor). The first differential pressure signal can represent the first differential pressure of the multiphase fluid 101 measured by the pressure sensors 104 a, 104 b (or first differential pressure sensor). The operations performed by the processor 705 can include receiving a second differential pressure signal from the second pair of pressure sensors 104 c, 104 d (or second differential pressure sensor). The second differential pressure signal can represent the second differential pressure of the multiphase fluid 101 measured by the pressure sensors 104 c, 104 d (or second differential pressure sensor). The operations performed by the processor 705 can include receiving a third differential pressure signal from the third pair of pressure sensors 104 e, 104 f (or third differential pressure sensor). The third differential pressure signal can represent the third differential pressure of the multiphase fluid 101 measured by the pressure sensors 104 e, 104 f (or third differential pressure sensor). The operations performed by the processor 705 can include receiving a fourth differential pressure signal from the fourth pair of pressure sensors 104 g, 104 h (or fourth differential pressure sensor). The fourth differential pressure signal can represent the fourth differential pressure of the multiphase fluid 101 measured by the pressure sensors 104 g, 104 h (or fourth differential pressure sensor). The operations performed by the processor 705 can include determining an overall density of the multiphase fluid 101 at least based on the specified height (h) and a difference between the third differential pressure and the fourth differential pressure. The operations performed by the processor 705 can include determining an overall flow rate of the multiphase fluid 101 at least based on the first differential pressure and the second differential pressure of the multiphase fluid 101 and the overall density of the multiphase fluid 101.

The mixture density of the multiphase fluid 101 can be determined in various ways. As a first example, the mixture density (ρ) of the multiphase fluid 101 can be determined by Equation 2d tailored for the apparatus 100H, where ΔP₁ is the third differential pressure measured by the third pair of pressure sensors 104 e, 104 f (or third differential pressure sensor), ΔP₂ is the fourth differential pressure measured by the fourth pair of pressure sensors 104 g, 104 h (or fourth differential pressure sensor), and h₂ is the specified height (h). As a second example, the mixture density (ρ) of the multiphase fluid 101 can be determined by Equation 5d tailored for the apparatus 100H, where ΔP₁ is the third differential pressure measured by the third pair of pressure sensors 104 e, 104 f (or third differential pressure sensor) and ΔP₂ is the fourth differential pressure measured by the fourth pair of pressure sensors 104 g, 104 h (or fourth differential pressure sensor).

The total flow rate of the multiphase fluid 101 can be determined in various ways. As a first example, the total flow rate of the multiphase fluid 101 can be determined using the first differential pressure by Equation 4a. As a second example, the total flow rate of the multiphase fluid 101 can be determined using the second differential pressure by Equation 4b. The operations performed by the processor 705 can include comparing the overall flow rates (m_(T1), m_(T2)) determined by Equations 4a and 4bc to determine an adjustment of the discharge coefficients (C_(d1), C_(d2))

In some implementations, a distance between the first pressure sensor 104 a and the second pressure sensor 104 b is substantially the same as a distance between the third pressure sensor 104 c and the fourth pressure sensor 104 d. In a specific example, the apparatus 100D has the following dimensions. The inner diameter (d) of the U-bend 102 (the first conduit 102 a, the second conduit 102 b, and the connecting conduit 102 c) is about 2.067 inches (NPS 2″ Schedule 40). The inner diameter (D) of the first pipe 150 a and the second pipe 150 b is about 3.068″ (NPS 3″ Schedule 40). The dimension A is the length of straight pipe (that is, pipe that is free of any bends) of the first pipe 150 a upstream of the first pressure sensor 104 a. The dimension A is also the length of straight pipe of the second pipe 150 b downstream of the fourth pressure sensor 104 d. The dimension A is a minimum length of straight pipe that ensures minimal interference of turbulence and friction which may affect the measurements taken by the first pair of pressure sensors 104 a, 104 b and the measurements taken by the second pair of pressure sensors 104 c, 104 d. In some implementations, the dimension A is at least five times the inner diameter (D) of the first pipe 150 a (A≥5D). The dimension B is the length of straight pipe of the first pipe 150 a downstream of the first pressure sensor 104 a and upstream of the first angled bend 153 a. The dimension B is also the length of straight pipe of the second pipe 150 b downstream of the second angled bend 153 b and upstream of the fourth pressure sensor 104 d. Similarly, the dimension B is a minimum length of straight pipe that ensures minimal interference of turbulence and friction which may affect the measurements taken by the first pair of pressure sensors 104 a, 104 b and the measurements taken by the second pair of pressure sensors 104 c, 104 d. In some implementations, the dimension B is at least five times the inner diameter (D) of the first pipe 150 a (B≥5D). The dimension C is the length of straight pipe of the first conduit 102 a downstream of the first angled bend 153 a and upstream of the second pressure sensor 104 b. The dimension C is also the length of straight pipe of the second conduit 102 b downstream of the third pressure sensor 104 c and upstream of the second angled bend 153 b. Similarly, the dimension C is a minimum length of straight pipe that ensures minimal interference of turbulence and friction which may affect the measurements taken by the first pair of pressure sensors 104 a, 104 b and the measurements taken by the second pair of pressure sensors 104 c, 104 d. In some implementations, the dimension C is at least five times the inner diameter (d) of the U-bend 102 (C≥5d). In some implementations, there is a length of straight pipe between the second pressure sensor 104 b and the fifth pressure sensor 104 e of at least two times the inner diameter (d) of the U-bend 102 (≥2d) along the longitudinal length of the first conduit 102 a of the U-bend 102. In some implementations, there is a length of straight pipe between the eighth pressure sensor 104 h and the third pressure sensor 104 c of at least two times the inner diameter (d) of the U-bend 102 (≥2d) along the longitudinal length of the second conduit 102 b of the U-bend 102. The dimension E is the length of straight pipe of the first conduit 102 a downstream of the sixth pressure sensor 104 f and upstream of the first curved bend 103 a. The dimension E is also the length of straight pipe of the second conduit 102 b downstream of the second curved bend 103 b and upstream of the seventh pressure sensor 104 g. Similarly, the dimension E is a minimum length of straight pipe that ensures minimal interference of turbulence and friction which may affect the measurements taken by the third pair of pressure sensors 104 e, 104 f and the measurements taken by the fourth pair of pressure sensors 104 g, 104 h. In some implementations, the dimension E is at least five times the inner diameter (d) of the U-bend 102 (E≥5d).

FIG. 1I depicts an apparatus 100I which can be used to measure a density, a flow rate, and a water-liquid ratio of a multiphase fluid 101. The apparatus 100I can be substantially similar to the apparatus 100A shown in FIG. 1A, the apparatus 100E shown in FIG. 1E, the apparatus 100G shown in FIG. 1G, and the apparatus 100H shown in FIG. 1H. For example, the apparatus 100I can be a combination of the apparatus 100G and the apparatus 100H. For example, the apparatus 100I can be a copy of the apparatus 100H that also includes dynamic pressure sensors 106 a, 106 b, 106 c (similar to the apparatus 100G) that are coupled to the U-bend 102. Thus, the apparatus 100I can perform all of the functions that can be performed by the apparatus 100G and also by the apparatus 100H.

Although shown in FIG. 1I as including three dynamic pressure sensors 106 a, 106 b, 106 c, the apparatus 100I can include fewer dynamic pressure sensors (such as two dynamic pressure sensors) or more dynamic pressure sensors (such as four dynamic pressure sensors or more than four dynamic pressure sensors). The dynamic pressure sensors 106 a, 106 b, 106 c measure a dynamic pressure of the multiphase fluid 101 flowing through the U-bend 102, as opposed to a static pressure of the multiphase fluid 101. The dynamic pressure sensors 106 a, 106 b, 106 c can be located at any locations along the U-bend 102. In some implementations, the dynamic pressure sensors 106 a, 106 b, 106 c are coupled to the first conduit 102 a of the U-bend 102. In some implementations, the dynamic pressure sensors 106 a, 106 b, 106 c are coupled to the second conduit 102 b of the U-bend 102. In some implementations, a first portion of the dynamic pressure sensors 106 a, 106 b, 106 c (for example, one or two of the dynamic pressure sensors 106 a, 106 b, 106 c) is coupled to the first conduit 102 a of the U-bend 102, and a remaining portion of the dynamic pressure sensors 106 a, 106 b, 106 c is coupled to the second conduit 102 b of the U-bend 102. The dynamic pressure sensors 106 a, 106 b, 106 c can be evenly or randomly distributed along the U-bend 102.

The dynamic pressure measurements taken by the dynamic pressure sensors 106 a, 106 b, 106 c can be cross-correlated (for example, by the processor 705) to determine a speed of sound of the multiphase fluid 101 flowing through the U-bend 102. The speed of sound of the multiphase fluid 101 and the overall density (ρ) of the multiphase fluid 101 can then be used to determine the water-liquid ratio of the multiphase fluid 101. For example, the speed of sound of the multiphase fluid 101 and the overall density (ρ) of the multiphase fluid 101 can be compared to a plotted graph of various phase curves (gas/oil, gas/aqueous, oil/aqueous) to determine the water-liquid ratio.

FIG. 1J depicts an apparatus 100J which can be used to measure a density, a flow rate, and a water-liquid ratio of a multiphase fluid 101. The apparatus 100J can be substantially similar to the apparatus 100H shown in FIG. 1H. In comparison, the apparatus 100J omits the connecting conduit (102 c shown in FIG. 1H). The apparatus 100J includes the U-bend 102, which includes the first conduit 102 a and the second conduit 102 b. The first conduit 102 a is connected to the second conduit 102 b by a curved bend 103. In some implementations, the first conduit 102 a, the second conduit 102 b, and the curved bend 103 are integrated, such that the U-bend 102 is a singular, unitary body, as opposed to parts that are disjointed and coupled together to form the U-bend 102. The curved bend 103 can be made of the same material as the first conduit 102 a and the second conduit 102 b, such that the friction experienced by the multiphase fluid 101 flowing through each of the components of the U-bend 102 is substantially the same.

As shown in FIG. 1J, the U-bend 102 can be connected to a flow pipe 150 that flows the multiphase fluid 101. The flow pipe 150 includes the first pipe 150 a and the second pipe 150 b. In some implementations, the U-bend 102 is connected to the flow pipe 150 by curved bends 153 c, 153 d. For example, the first conduit 102 a of the U-bend 102 is connected to the first pipe 150 a of the flow pipe 150 by the first curved bend 153 c. For example, the second conduit 102 b of the U-bend 102 is connected to the second pipe 150 b of the flow pipe 150 by the second curved bend 153 d. In some implementations, the first conduit 102 a, the second conduit 102 b, the curved bend 103, the first curved bend 153 c, the second curved bend 153 d, and the flow pipe 150 are integrated, such that the U-bend 102 and the flow pipe 150 form a singular, unitary body, as opposed to parts that are disjointed and coupled together.

The apparatus 100J can include a first differential pressure sensor 108 a coupled to the first pipe 150 a and to the first conduit 102 a via pressure ports 105 a and 105 b, respectively. The first differential pressure sensor 108 a is configured to measure a first differential pressure of the multiphase fluid 101 flowing across the first curved bend 153 c. The pressure ports 105 a and 105 b are separated by the specified differential pressure height, h, with respect to gravity. The apparatus 100B can include a second differential pressure sensor 108 b coupled to the second conduit 102 b and to the second pipe 150 b via pressure ports 105 c and 105 d, respectively. The second differential pressure sensor 108 b is configured to measure a second differential pressure of the multiphase fluid 101 flowing across the second curved bend 153 d. The pressure ports 105 c and 105 d are separated by the specified differential pressure height, h, with respect to gravity. As shown in FIG. 1F, the pressure port 105 a does not necessarily need to be at the same vertical height as pressure port 105 d, and the pressure port 105 b does not necessarily need to be at the same vertical height as pressure port 105 c. Instead, the vertical distance between pressure port 105 a and pressure port 105 b should match the vertical distance between pressure port 105 c and pressure port 105 d. The height difference between the pressure ports 105 a, 105 b of the first differential pressure sensor 108 a and the height difference between the pressure ports 105 c, 105 d of the second differential pressure sensor 108 b are the same, and the horizontal distance between the pressure ports 105 a, 105 b of the first differential pressure sensor 108 a and the horizontal distance between the pressure ports 105 c, 105 d of the second differential pressure sensor 108 b are the same, such that the frictional components of the pressure drops experienced by the multiphase fluid 101 flowing across the first curved bend 153 c (measured by the first differential pressure sensor 108 a) and through the second curved bend 153 d (measured by the second differential pressure sensor 108 b) are substantially the same.

The apparatus 100J can include a third differential pressure sensor 108 c coupled to the first conduit 102 a via pressure ports 105 e and 105 f. The third differential pressure sensor 108 c is configured to measure a third differential pressure of the multiphase fluid 101 flowing through the U-bend 102. The pressure ports 105 e and 105 f are separated by the specified differential pressure height, h, with respect to gravity. The apparatus 100J can include a fourth differential pressure sensor 108 d coupled to the second conduit 102 b via pressure ports 105 g and 105 h. The fourth differential pressure sensor 108 d is configured to measure a fourth differential pressure of the multiphase fluid 101 flowing through the U-bend 102. The pressure ports 105 g and 105 h are separated by the specified differential pressure height, h, with respect to gravity. As shown in FIG. 1J, the pressure port 105 e does not necessarily need to be at the same vertical height as pressure port 105 h, and the pressure port 105 f does not necessarily need to be at the same vertical height as pressure port 105 g. Instead, the vertical distance between pressure port 105 e and pressure port 105 f should match the vertical distance between pressure port 105 g and pressure port 105 h. The height difference between the pressure ports 105 e, 105 f of the third differential pressure sensor 108 c and the height difference between the pressure ports 105 g, 105 h of the fourth differential pressure sensor 108 d are the same, such that the frictional components of the pressure drops experienced by the multiphase fluid 101 flowing through the first conduit 102 a (measured by the third differential pressure sensor 108 c) and through the second conduit 102 b (measured by the fourth differential pressure sensor 108 d) are substantially the same.

Although not shown in FIG. 1J, the first, second, third, and fourth differential pressure sensors 108 a, 108 b, 108 c, 108 d can be coupled to a computer (such as the computer 700), to determine characteristics (such as a mixture density, a total mass flow rate, an aqueous phase volume fraction, a gas phase volume fraction, or any combination of these) of the multiphase fluid 101 flowing through the U-bend 102 (for example, using Equations 1, 1a, 1b, 2a, 2b, 2c, 2d, 2e, 2f, 2g, 2h, 3a, 3b, 4a, 4b, 5a, 5b, 5c, 5d, or any combination of these).

Any of the apparatuses 100A, 100B, 100C, 100D, 100E, 100F, 100G, 100H, 100I, or 100J can be combined with any of the following components to obtain additional information about the multiphase fluid 101: a flowmeter (such as a Venturi-type meter), a water-cut meter that measures water fraction or water-liquid ratio, a dynamic pressure sensor, a temperature sensor, or a strain gauge. In cases where a dynamic pressure sensor is included, the dynamic pressure sensor can be located on the same pressure tap as any of the pressure sensors already included in the apparatus (100A, 100B, 100C, 100D, 100E, 100F, 100G, 100H, 100I, or 100J).

FIG. 2A is a schematic diagram depicting various example pressure tap arrangements. Each of the views (i), (ii), and (iii) are cross-sectional views of a pipe (for example, the first conduit 102 a or the second conduit 102 b). Any of the pressure sensors (104 a, 104 b, 104 c, 104 d, 104 e, 104 f, 104 a′, 104 b′, 104 c′, 104 d′, 104 e′, 104 f′) can have the form shown in view (i), (ii), or (iii). In view (i), the pressure sensor includes a single pressure tap. In view (ii), the pressure sensor includes various pressure taps (for example, four) that are distributed around a circumference of the conduit. In view (ii), the pressure sensor can obtain four simultaneous readings (one at each of the four pressure taps) and average the four simultaneous readings to obtain an averaged pressure reading for that particular axial location along the conduit. Although shown in view (ii) as having four pressure taps, the pressure sensors can have fewer pressure taps (such as two or three) or more pressure taps (such as five or more than five). In view (iii), the pressure sensor includes a circumferential pressure tap that fully spans the entire circumference of the conduit. The circumferential pressure tap shown in view (iii) may be the most accurate out of all the views (i), (ii), and (iii) shown in FIG. 2A.

FIG. 2B is a schematic diagram depicting various example piping arrangements. Each of the views (i), (ii), (iii), (iv), and (v) depict piping arrangements that can used in addition to or alternatively to the piping arrangements already described with respect to the apparatuses 100A, 100B, 100C, 100D, 100E, 100F, 100G, 100H, 100I, and 100J. ΔP₁ is a first differential pressure of the multiphase fluid 101 flowing through the respective piping arrangement. ΔP₂ is a second differential pressure of the multiphase fluid 101 flowing through the respective piping arrangement. ΔP₁ can be measured by a pair of pressure sensors or a single differential pressure sensor. Similarly, ΔP₂ can be measured by a pair of pressure sensors or a single differential pressure sensor. The first and second differential pressures (ΔP₁, ΔP₂) measured from any of these piping arrangements can be used, for example, to determine characteristics (such as a mixture density, a total mass flow rate, an aqueous phase volume fraction, a gas phase volume fraction, or any combination of these) of the multiphase fluid 101 flowing through the respective piping arrangement (for example, using Equations 1, 1a, 1b, 2a, 2b, 2c, 2d, 2e, 2f, 2g, 2h, 3a, 3b, or any combination of these).

As an example, view (i) can be considered a base case that generally conforms to the piping arrangements of apparatuses 100A and 100B shown in FIGS. 1A and 1B, respectively. The piping arrangement of view (ii) is substantially similar to the piping arrangement of view (i), but the first and second conduits 102 a, 102 b of view (ii) are not necessarily straight pipes. In view (iii), the connecting conduit (102 c) is omitted, and the first and second conduits 102 a, 102 b are connected by an angled bend. The piping arrangement of view (iv) is substantially similar to the piping arrangement of view (i), but is flipped vertically, such that the multiphase fluid 101 flowing through the piping arrangement of view (iv) flows downward before flowing upward. The piping arrangement of view (v) is substantially similar to the piping arrangement of view (i), but the first and second conduits 102 a, 102 b are coiled (that is, includes coiled piping).

FIG. 2C is a schematic diagram depicting various example piping arrangements. Each of the views (i), (ii), (iii), (iv), and (v) depict piping arrangements that can used in addition to or alternatively to the piping arrangements already described with respect to the apparatuses 100A, 100B, 100C, 100D, 100E, 100F, 100G, 100H, 100I, and 100J. ΔP₁ is a first differential pressure of the multiphase fluid 101 flowing through the respective piping arrangement across a first bend (curved or angled). ΔP₂ is a second differential pressure of the multiphase fluid 101 flowing through the respective piping arrangement across a second bend (curved or angled). ΔP₁ can be measured by a pair of pressure sensors or a single differential pressure sensor. Similarly, ΔP₂ can be measured by a pair of pressure sensors or a single differential pressure sensor. The first and second differential pressures (ΔP₁, ΔP₂) measured from any of these piping arrangements can be used, for example, to determine characteristics (such as a mixture density, a total mass flow rate, or both) of the multiphase fluid 101 flowing through the respective piping arrangement (for example, using Equations 4a, 4b, 5a, 5b, 5c, 5d, or any combination of these).

As an example, view (i) can be considered a base case that generally conforms to the piping arrangements of apparatuses 100E and 100F shown in FIGS. 1E and 1F, respectively. The piping arrangement of view (ii) is substantially similar to the piping arrangement of view (i), but the first and second conduits 102 a, 102 b of view (ii) deviate from vertically oriented pipes (that is, they are disposed at non-zero angles with respect to a vertical) and are not necessarily straight pipes. Similarly, in view (iii), the first and second conduits 102 a, 102 b deviate from vertically oriented pipes, but are straight pipes that are connected with angled bends. The piping arrangement of view (iv) is substantially similar to the piping arrangement of view (i), but is flipped vertically, such that the multiphase fluid 101 flowing through the piping arrangement of view (iv) flows downward before flowing upward. The piping arrangement of view (v) is substantially similar to the piping arrangement of view (i), but the first and second conduits 102 a, 102 b are coiled (that is, includes coiled piping).

FIG. 2D is a schematic diagram depicting various example piping arrangements. Each of the views (i), (ii), (iii), and (iv) depict piping arrangements that can used in addition to or alternatively to the piping arrangements already described with respect to the apparatuses 100A, 100B, 100C, 100D, 100E, 100F, 100G, 100H, 100I, and 100J. ΔP₁ is a first differential pressure of the multiphase fluid 101 flowing through the respective piping arrangement across a first bend (curved or angled). ΔP₂ is a second differential pressure of the multiphase fluid 101 flowing through the respective piping arrangement across a second bend (curved or angled). ΔP₁ can be measured by a pair of pressure sensors or a single differential pressure sensor. Similarly, ΔP₂ can be measured by a pair of pressure sensors or a single differential pressure sensor. ΔP₃ is a third differential pressure of the multiphase fluid 101 flowing through the first conduit 102 a of the respective piping arrangement. ΔP₄ is a fourth differential pressure of the multiphase fluid 101 flowing through the second conduit 102 b of the respective piping arrangement. ΔP₃ can be measured by a pair of pressure sensors or a single differential pressure sensor. Similarly, ΔP₄ can be measured by a pair of pressure sensors or a single differential pressure sensor. The first, second, third, and fourth differential pressures (ΔP₁, ΔP₂, ΔP₃, ΔP₄) measured from any of these piping arrangements can be used, for example, to determine characteristics (such as a mixture density, a total mass flow rate, an aqueous phase volume fraction, a gas phase volume fraction, or any combination of these) of the multiphase fluid 101 flowing through the respective piping arrangement (for example, using Equations 1, 1a, 1b, 2a, 2b, 2c, 2d, 2e, 2f, 2g, 2h, 3a, 3b, 4a, 4b, 5b, 5c, 5d, or any combination of these).

As an example, view (i) can be considered a base case that generally conforms to the piping arrangements of apparatuses 100H and 100J shown in FIGS. 1H and 1J, respectively. The piping arrangement of view (ii) is substantially similar to the piping arrangement of view (i), but the first and second conduits 102 a, 102 b of view (ii) deviate from vertically oriented pipes and are connected with angled bends. In the piping arrangement of view (iii), the first and second pipes 150 a, 150 b of the flow pipe 150 are not in-line with one another. The piping arrangement of view (iv) is substantially similar to the piping arrangement of view (i), but is flipped vertically, such that the multiphase fluid 101 flowing through the piping arrangement of view (iv) flows downward before flowing upward.

FIG. 2E is a schematic diagram depicting various example piping arrangements. Each of the views (i), (ii), (iii), and (iv) depict piping arrangements that can used in addition to or alternatively to the piping arrangements already described with respect to the apparatuses 100A, 100B, 100C, 100D, 100E, 100F, 100G, 100H, 100I, and 100J. ΔP₁ is a first differential pressure of the multiphase fluid 101 flowing through the respective piping arrangement across a first bend (curved or angled). ΔP₂ is a second differential pressure of the multiphase fluid 101 flowing through the respective piping arrangement across a second bend (curved or angled). ΔP₁ can be measured by a pair of pressure sensors or a single differential pressure sensor. Similarly, ΔP₂ can be measured by a pair of pressure sensors or a single differential pressure sensor. ΔP₃ is a third differential pressure of the multiphase fluid 101 flowing through the first conduit 102 a of the respective piping arrangement. ΔP₄ is a fourth differential pressure of the multiphase fluid 101 flowing through the second conduit 102 b of the respective piping arrangement. ΔP₃ can be measured by a pair of pressure sensors or a single differential pressure sensor. Similarly, ΔP₄ can be measured by a pair of pressure sensors or a single differential pressure sensor. The first, second, third, and fourth differential pressures (ΔP₁, ΔP₂, ΔP₃, ΔP₄) measured from any of these piping arrangements can be used, for example, to determine characteristics (such as a mixture density, a total mass flow rate, an aqueous phase volume fraction, a gas phase volume fraction, or any combination of these) of the multiphase fluid 101 flowing through the respective piping arrangement (for example, using Equations 1, 1a, 1b, 2a, 2b, 2c, 2d, 2e, 2f, 2g, 2h, 3a, 3b, 4a, 4b, 5b, 5c, 5d, or any combination of these).

The piping arrangement of view (i) is substantially similar to the piping arrangements of apparatuses 100H and 100J shown in FIGS. 1H and 1J, respectively, but the first and second conduits 102 a, 102 b of view (i) are not necessarily straight pipes. In the piping arrangement of view (ii), the first and second conduits 102 a, 102 b form a loop. The piping arrangement of view (iii) is substantially similar to the piping arrangement of view (i) of FIG. 2D, but the first and second conduits 102 a, 102 b are coiled (that is, includes coiled piping). The piping arrangement of view (iv) is substantially similar to the piping arrangement of view (ii), but the first and second conduits 102 a, 102 b form an angled loop with straight portions of piping in between angled bends. Further, in the piping arrangement of view (iv), the first and second pipes 150 a, 150 b of the flow pipe 150 are not in-line with one another

FIG. 2F is a schematic diagram depicting various example sensor arrangements. Any of the sensor arrangements shown in views (i), (ii), (iii), (iv), and (v) can be applied (either alone or in combination) in addition to or alternatively to the sensor arrangements already described with respect to the apparatuses 100A, 100B, 100C, 100D, 100E, 100F, 100G, 100H, 100I, and 100J.

As an example, view (i) can be considered a base case that generally conforms to the differential pressure sensor arrangements of apparatuses 100B, 100F, and 100J shown in FIGS. 1B, 1F, and 1J, respectively. The pressure sensor arrangement of view (ii) includes two pressure sensors which measure static pressure at their respective locations. A difference between the measured static pressures can be determined as a differential pressure (similar to the differential pressure sensor). The pressure sensor arrangement of view (ii) is substantially similar to the pressure sensor arrangements shown in the apparatuses 100A, 100C, 100D, 100E, 100G, 100H, and 100I shown in FIGS. 1A, 1C, 1D, 1E, 1G, 1H, and 1I, respectively. The differential pressure sensor arrangement of view (iii) is substantially similar to the view (i) but includes redundant differential pressure sensors in a parallel configuration across the same length of piping, which is substantially similar to the differential pressure sensor arrangement shown in apparatus 100C of FIG. 1C. The differential pressure sensor arrangement of view (iv) is substantially similar to the differential pressure sensor arrangement of view (i), but includes additional differential pressure sensors in a series configuration across a length of piping, which is substantially similar to the differential pressure sensor arrangement shown in apparatus 100D of FIG. 1D. The sensor arrangement of view (v) includes multiple strain sensors (such as strain gauges) disposed along a length of piping. The measured strain can be correlated to pressure, such that similar calculations can be completed to determine characteristics of the flowing multiphase fluid 101.

FIG. 3A is a flowchart depicting an example multiphase fluid flow model 300A. Any one of the apparatuses 100H, 100I, or 100J can, for example, be used to implement the model 300A to characterize a multiphase fluid, such as the multiphase fluid 101. At block 301 a, pressure and temperature of the multiphase fluid 101 is determined (for example, using pressure and temperature sensors). At block 303 a, various properties (such as the density, viscosity, and speed of sound) of the phases (gas phase 101 a, aqueous phase 101 b′, and oil phase 101 b″) of the multiphase fluid 101 are determined. For example, the gas phase density (ρ_(g)) of the gas phase 101 a, the aqueous phase density (ρ_(a)) of the aqueous phase 101 b′, and the oil phase density (ρ_(o)) of the oil phase 101 b″ are determined at block 303 a. At block 305 a, a total flow rate (m_(T)) of the multiphase fluid 101 is determined (for example, using differential pressure sensors as in the apparatus 100D). At block 307 a, a mixture density (ρ) of the multiphase fluid 101 is determined (for example, using differential pressure sensors as in the apparatus 100A). At block 309 a, a water-liquid ratio (R_(WL)) of the multiphase fluid 101 is determined (for example, using dynamic pressure sensors as in the apparatus 100E)). At block 311 a, a gas volume fraction (α_(g)) is determined at least based on the densities of the various phases (determined at block 303 a), the mixture density (ρ) of the multiphase fluid 101 (determined at block 307 a), and the water-liquid ratio (R_(WL)) of the multiphase fluid 101 (determined at block 309 a): ρ=(ρ_(a)×R_(WL)+ρ_(o)×(1−R_(WL)))×(1−α_(g))+ρ_(g)×α_(g). At block 313 a, a gas flow rate (Q_(g)) of the gas phase 101 a is determined based on the gas volume fraction (α_(g)) and the total volume flow rate (Q) of the multiphase fluid 101: Q_(g)=α_(g)×Q. At block 313 a, an aqueous phase flow rate (Q_(a)) of the aqueous phase 101 b′ is determined based on the water-liquid ratio (R_(WL)), the gas volume fraction (α_(g)), and the overall flow rate (Q) of the multiphase fluid 101: Q_(a)=R_(WL)×(1−α_(g))×Q. At block 313 a, an oil phase flow rate (Q_(o)) of the oil phase 101 b″ is determined based on the water-liquid ratio (R_(WL)), the gas volume fraction (α_(g)), and the total volume flow rate (Q) of the multiphase fluid 101: Q_(o)=(1−R_(WL))×(1−α_(g))×Q.

FIG. 3B is a flowchart depicting an example multiphase fluid flow model 300B. Any one of the apparatuses 100A, 100B, 100C, 100D, 100H, 100I, or 100J can, for example, be used to implement the model 300B to characterize a multiphase fluid, such as the multiphase fluid 101. At block 301 b, a first differential pressure (ΔP₁) and a second differential pressure (ΔP₂) of the flowing multiphase fluid 101 is measured. The first differential pressure (ΔP₁) is measured in a first direction across a specified differential pressure height (h). The second differential pressure (ΔP₂) is measured in a second direction (different from the first direction) across the specified differential pressure height (h). The first and second directions are such that the multiphase fluid 101 flows generally upward across the specified differential pressure height (h) in one of the first or second directions and flows generally downward across the specified differential pressure height (h) in the other direction. For example, the multiphase fluid 101 flows generally upward across the specified differential pressure height (h) in the first direction and flows generally downward across the specified differential pressure height (h) in the second direction. For example, the multiphase fluid 101 flows generally downward across the specified differential pressure height (h) in the first direction and flows generally upward across the specified differential pressure height (h) in the second direction.

At block 303 b, a frictional pressure loss (ΔP_(f)) of the flowing multiphase fluid 101 is calculated based on the first and second differential pressures (ΔP₁, ΔP₂) measured at block 301 b. The frictional pressure loss (ΔP_(f)) can be calculated at block 303 b, for example, by Equation 2h. At block 305 b, a mixture density (ρ) of the multiphase fluid 101 is calculated based on the first and second differential pressures (ΔP₁, ΔP₂) measured at block 301 b. The mixture density (ρ) can be calculated at block 305 b, for example, by Equation 2d. At block 307 b, a gas phase volume fraction (α_(g)) of the multiphase fluid 101 is calculated based on the mixture density (ρ) determined at block 305 b. The gas phase volume fraction (α_(g)) can be calculated at block 307 b, for example, by Equation 3b. In some implementations, gas phase volume fraction (α_(g)) can be calculated at block 307 b based on other known or measured parameters, such as density of the liquid phase (ρ_(l)) and density of the gas phase (ρ_(g)). In such implementations, the gas phase volume fraction (α_(g)) can be calculated by Equation 6a.

$\begin{matrix} {\alpha_{g} = \frac{\rho - \rho_{l}}{\rho_{g} - \rho_{l}}} & \left( {6a} \right) \end{matrix}$

In some implementations, a total mass flow rate (m_(T)) of the multiphase fluid 101 is determined. The total mass flow rate (m_(T)) can be determined, for example, by a mass flowmeter and/or by any of Equations 3a, 4a, or 4b. The multiphase fluid flow model 300B assumes that a slip ratio (S) of the multiphase fluid 101 is substantially equal to 1. The slip ratio (S) is the ratio of the velocity of the gas phase to the velocity of the liquid phase (including the aqueous phase and the oil phase) in the flowing multiphase fluid 101.

FIG. 3C is a flowchart depicting an example multiphase fluid flow model 300C. Any one of the apparatuses 100A, 100B, 100C, 100D, 100H, 100I, or 100J can, for example, be used to implement the model 300C to characterize a multiphase fluid, such as the multiphase fluid 101. At block 301 c, a first differential pressure (ΔP₁) and a second differential pressure (ΔP₂) of the flowing multiphase fluid 101 is measured. The first differential pressure (ΔP₁) is measured in a first direction across a specified differential pressure height (h). The second differential pressure (ΔP₂) is measured in a second direction (different from the first direction) across the specified differential pressure height (h). The first and second directions are such that the multiphase fluid 101 flows generally upward across the specified differential pressure height (h) in one of the first or second directions and flows generally downward across the specified differential pressure height (h) in the other direction. For example, the multiphase fluid 101 flows generally upward across the specified differential pressure height (h) in the first direction and flows generally downward across the specified differential pressure height (h) in the second direction. For example, the multiphase fluid 101 flows generally downward across the specified differential pressure height (h) in the first direction and flows generally upward across the specified differential pressure height (h) in the second direction.

At block 303 c, a frictional pressure loss (ΔP_(f)) of the flowing multiphase fluid 101 is calculated based on the first and second differential pressures (ΔP₁, ΔP₂) measured at block 301 c. The frictional pressure loss (ΔP_(f)) can be calculated at block 303 c, for example, by Equation 2h. At block 305 c, an apparent mixture density (ρ_(app)) of the multiphase fluid 101 is estimated based on the first and second differential pressures (ΔP₁, ΔP₂) measured at block 301 c. The apparent mixture density (ρ_(app)) can be estimated at block 305 c, for example, by Equation 2d. At block 307 c, a void fraction (λ_(g)) of the multiphase fluid 101 is calculated based on the apparent mixture density (ρ_(app)) estimated at block 305 c. The void fraction (λ_(g)) can be calculated at block 307 c, for example, by Equation 6b.

$\begin{matrix} {\lambda_{g} = \frac{\rho_{app} - \rho_{l}}{\rho_{g} - \rho_{l}}} & \left( {6b} \right) \end{matrix}$

At block 309 c, various operating conditions and/or properties (such as pressure, temperature, and salinity) about the flowing multiphase fluid 101 are determined. The conditions and/or properties determined at block 309 c can be measured, for example, by external sensors (such as pressure sensors, temperature sensors, and brine samplers). At block 311 c, a water-liquid ratio (R_(WL)) of the flowing multiphase fluid 101 is determined. The water-liquid ratio (R_(WL)) can be determined at block 311 c can be measured, for example, by external sensors (such as dynamic pressure sensors, as in apparatus 100G shown in FIG. 1G). At block 313 c, densities of the various phases (such as density of the gas phase (ρ_(g)), density of the aqueous phase (ρ_(a)), and density of the oil phase (ρ_(o))) of the flowing multiphase fluid 101 are determined. The densities of the various phases of the flowing multiphase fluid 101 can be determined at block 313 c, for example, based on the operating conditions (pressure and temperature) of the flowing multiphase fluid 101 determined at block 309 c. At block 313 c, a density of the liquid phase (ρ_(l)), including the aqueous phase and the oil phase) is determined. The density of the liquid phase (ρ_(l)) can be determined at block 313 c based on the densities of the various phases of the flowing multiphase fluid 101 and the water-liquid ratio (R_(WL)) determined at block 311 c. The density of the liquid phase (ρ_(l)) can be determined at block 313 c by Equation 6c.

ρ_(l)=ρ_(a) R _(WL)+ρ_(o)(1−R _(WL))  (6c)

In some implementations, the void fraction (λ_(g)) can be calculated at block 307 c based on the density of the liquid phase (ρ_(l)) and the density of the gas phase (ρ_(g)) determined at block 313 c.

At block 315 c, a gas phase volume fraction (α_(g)) of the multiphase fluid 101 is calculated based on the void fraction (λ_(g)) determined at block 307 c and slip ratio (S) of the flowing multiphase fluid 101. The slip ratio (S) used at block 315 c to determine the gas phase volume fraction (α_(g)) can, for example, be measured by an external sensor or assumed to be substantially equal to 1 as an initial guess. At block 317 c, a mixture density (ρ) of the multiphase fluid 101 is calculated based on the gas phase volume fraction (α_(g)) determined at block 315 c, the density of the gas phase (ρ_(g)) determined at block 313 c, and the density of the liquid phase (ρ_(l)) determined at block 313 c. The mixture density (ρ) can be determined at block 317 c by Equation 6d.

ρ=α_(g)ρ_(g)+(1−α_(g))ρ_(l)  (6d)

In some implementations, a total mass flow rate (m_(T)) of the multiphase fluid 101 is determined. The total mass flow rate (m_(T)) can be determined, for example, by a mass flowmeter and/or by any of Equations 3a, 4a, or 4b. In some implementations, the blocks in the dotted box of the multiphase fluid flow model 300C (blocks 315 c, 317 c, determining the total mass flow rate (m_(T)), and determining the slip ratio (S)) are performed simultaneously (for example, in parallel). In some implementations, the blocks in the dotted box of the multiphase fluid flow model 300C (blocks 315 c, 317 c, determining the total mass flow rate (m_(T)), and determining the slip ratio (S)) are performed iteratively (for example, in series) until the calculated values converge to final values.

FIG. 3D is a flowchart depicting an example multiphase fluid flow model 300D. Any one of the apparatuses 100A, 100B, 100C, 100D, 100H, 100I, or 100J can, for example, be used to implement the model 300D to characterize a multiphase fluid, such as the multiphase fluid 101. At block 301 d, a first differential pressure (ΔP₁) and a second differential pressure (ΔP₂) of the flowing multiphase fluid 101 is measured. The first differential pressure (ΔP₁) is measured in a first direction across a specified differential pressure height (h). The second differential pressure (ΔP₂) is measured in a second direction (different from the first direction) across the specified differential pressure height (h). The first and second directions are such that the multiphase fluid 101 flows generally upward across the specified differential pressure height (h) in one of the first or second directions and flows generally downward across the specified differential pressure height (h) in the other direction. For example, the multiphase fluid 101 flows generally upward across the specified differential pressure height (h) in the first direction and flows generally downward across the specified differential pressure height (h) in the second direction. For example, the multiphase fluid 101 flows generally downward across the specified differential pressure height (h) in the first direction and flows generally upward across the specified differential pressure height (h) in the second direction.

At block 303 d, a frictional pressure loss (ΔP_(f)) of the flowing multiphase fluid 101 is calculated based on the first and second differential pressures (ΔP₁, ΔP₂) measured at block 301 d. The frictional pressure loss (ΔP_(f)) can be calculated at block 303 d, for example, by Equation 2h. At block 305 d, an apparent mixture density (ρ_(app)) of the multiphase fluid 101 is estimated based on the first and second differential pressures (ΔP₁, ΔP₂) measured at block 301 d. The apparent mixture density (ρ_(app)) can be estimated at block 305 d, for example, by Equation 2d.

At block 307 d, various operating conditions and/or properties (such as pressure, temperature, and salinity) about the flowing multiphase fluid 101 are determined. The conditions and/or properties determined at block 307 d can be measured, for example, by external sensors (such as pressure sensors, temperature sensors, and brine samplers). At block 309 d, densities of the various phases (such as density of the gas phase (ρ_(g)), density of the aqueous phase (ρ_(a)), and density of the oil phase (ρ_(o))) of the flowing multiphase fluid 101 are determined. The densities of the various phases of the flowing multiphase fluid 101 can be determined at block 309 d, for example, based on the operating conditions (pressure and temperature) of the flowing multiphase fluid 101 determined at block 307 d. At block 311 d, a gas phase volume fraction (α_(g)) of the flowing multiphase fluid 101 is determined. The gas phase volume fraction (α_(g)) determined at block 311 d can be measured, for example, by an external sensor (such as an ultrasonic sensor).

At block 313 d, a void fraction (λ_(g)) of the multiphase fluid 101 is calculated based on the gas volume fraction fraction (α_(g)) determined at block 311 d and slip ratio (S) of the flowing multiphase fluid 101. The slip ratio (S) used at block 313 d to determine the void fraction (λ_(g)) can, for example, be measured by an external sensor or assumed to be substantially equal to 1 as an initial guess. At block 315 d, a density of the liquid phase (ρ_(l)), including the aqueous phase and the oil phase) is determined. The density of the liquid phase (ρ_(l)) can be determined at block 315 d based on the apparent mixture density (ρ_(app)) determined at block 305 d, the densities of the various phases of the flowing multiphase fluid 101 determined at block 309 d, and the void fraction (λ_(g)) determined at block 313 d. The density of the liquid phase (ρ_(l)) can be calculated at block 315 d, for example, by Equation 6e.

$\begin{matrix} {\rho_{l} = \frac{\rho_{app} - {\lambda_{g}\rho_{g}}}{1 - \lambda_{g}}} & \left( {6e} \right) \end{matrix}$

At block 317 d, a mixture density (ρ) of the multiphase fluid 101 is calculated based on the density of the gas phase (ρ_(g)) determined at block 309 d, the gas phase volume fraction (α_(g)) determined at block 311 d, and the density of the liquid phase (ρ_(l)) determined at block 315 d. The mixture density (ρ) can be determined at block 317 c by Equation 6d.

In some implementations, a total mass flow rate (m_(T)) of the multiphase fluid 101 is determined. The total mass flow rate (m_(T)) can be determined, for example, by a mass flowmeter and/or by any of Equations 3a, 4a, or 4b. In some implementations, the blocks in the dotted box of the multiphase fluid flow model 300D (blocks 313 d, 315 d, 317 d, determining the total mass flow rate (m_(T)), and determining the slip ratio (S)) are performed simultaneously (for example, in parallel). In some implementations, the blocks in the dotted box of the multiphase fluid flow model 300D (blocks 313 d, 315 d, 317 d, determining the total mass flow rate (m_(T)), and determining the slip ratio (S)) are performed iteratively (for example, in series) until the calculated values converge to final values. At block 319 d, a water-liquid ratio (R_(WL)) of the multiphase fluid 101 is determined. The water-liquid ratio (R_(WL)) can be determined at block 319 d by re-arranging Equation 6c to solve for

${R_{WL}\left( {R_{WL} = \frac{\rho_{l} - \rho_{o}}{\rho_{a} - \rho_{o}}} \right)}.$

FIG. 3E is a flowchart depicting an example multiphase fluid flow model 300E. Any one of the apparatuses 100E, 100F, 100H, 100I, or 100J can, for example, be used to implement the model 300E to characterize a multiphase fluid, such as the multiphase fluid 101. At block 301 e, a first differential pressure (ΔP₁) and a second differential pressure (ΔP₂) of the flowing multiphase fluid 101 is measured. The first differential pressure (ΔP₁) is measured across a first bend in piping. The second differential pressure (ΔP₂) is measured across a second bend in piping. The multiphase fluid 101 flows generally upward across either of the first or second bends and flows generally downward across the other of the first or second bends. For example, the multiphase fluid 101 flows generally upward across the first bend and flows generally downward across the second bend. For example, the multiphase fluid 101 flows generally downward across the first bend and flows generally upward across the second bend. Examples are provided with respect to apparatuses 100E, 100F, 100H, 100I, and 100J shown in FIGS. 1E, 1F, 1H, 1I, and 1J, respectively). At block 303 e, either one of the first differential pressure (ΔP₁) or the second differential pressure (ΔP₂) measured at block 301 e is selected for further calculations in the model 300E.

At block 305 e, a void fraction (λ_(g)) of the multiphase fluid 101 is determined. The void fraction (λ_(g)) determined at block 305 e can be measured for example, by an external sensor (such as an ultrasonic sensor). At block 307 e, various operating conditions and/or properties (such as pressure, temperature, and salinity) about the flowing multiphase fluid 101 are determined. The conditions and/or properties determined at block 307 e can be measured, for example, by external sensors (such as pressure sensors, temperature sensors, and brine samplers). At block 309 e, a water-liquid ratio (R_(WL)) of the flowing multiphase fluid 101 is determined. The water-liquid ratio (R_(WL)) can be determined at block 309 e can be measured, for example, by external sensors (such as dynamic pressure sensors, as in apparatus 100G shown in FIG. 1G). At block 311 e, densities of the various phases (such as density of the gas phase (ρ_(g)), density of the aqueous phase (ρ_(a)), and density of the oil phase (ρ_(o))) of the flowing multiphase fluid 101 are determined. The densities of the various phases of the flowing multiphase fluid 101 can be determined at block 311 e, for example, based on the operating conditions (pressure and temperature) of the flowing multiphase fluid 101 determined at block 307 e. At block 311 e, a density of the liquid phase (ρ_(l), including the aqueous phase and the oil phase) is determined. The density of the liquid phase (ρ_(l)) can be determined at block 311 e based on the densities of the various phases of the flowing multiphase fluid 101 and the water-liquid ratio (R_(WL)) determined at block 309 e. The density of the liquid phase (ρ_(l)) can be determined at block 311 e by Equation 6c. In some implementations, the void fraction (α_(g)) can be calculated at block 305 e based on the density of the liquid phase (ρ_(l)) and the density of the gas phase (ρ_(g)) determined at block 311 e.

At block 313 e, a gas phase volume fraction (α_(g)) of the multiphase fluid 101 is calculated based on the void fraction (λ_(g)) determined at block 305 e and slip ratio (S) of the flowing multiphase fluid 101. The slip ratio (S) used at block 313 e to determine the gas phase volume fraction (α_(g)) can, for example, be measured by an external sensor or assumed to be substantially equal to 1 as an initial guess. At block 315 e, a mixture density (ρ) of the multiphase fluid 101 is calculated based on the gas phase volume fraction (α_(g)) determined at block 313 e, the density of the gas phase (ρ_(g)) determined at block 311 e, and the density of the liquid phase (ρ_(l)) determined at block 311 e. The mixture density (ρ) can be determined at block 315 e by Equation 6d. At block 317 e, a total mass flow rate (m_(T)) of the multiphase fluid 101 is determined. The total mass flow rate (m_(T)) can be determined at block 317 e, for example, by a mass flowmeter and/or by any of Equations 4a or 4b based on the selected differential pressure at block 303 e. In some implementations, the blocks in the dotted box of the multiphase fluid flow model 300E (blocks 313 e, 315 e, 317 e, and determining the slip ratio (S)) are performed simultaneously (for example, in parallel). In some implementations, the blocks in the dotted box of the multiphase fluid flow model 300E (blocks 313 e, 315 e, 317 e, and determining the slip ratio (S)) are performed iteratively (for example, in series) until the calculated values converge to final values.

FIG. 3F is a flowchart depicting an example multiphase fluid flow model 300F. Any one of the apparatuses 100H, 100I, or 100J can, for example, be used to implement the model 300F to characterize a multiphase fluid, such as the multiphase fluid 101. At block 301 f, a first differential pressure (ΔP₁) and a second differential pressure (ΔP₂) of the flowing multiphase fluid 101 is measured. The first differential pressure (ΔP₁) is measured across a first bend in piping. The second differential pressure (ΔP₂) is measured across a second bend in piping. The multiphase fluid 101 flows generally upward across either of the first or second bends and flows generally downward across the other of the first or second bends. For example, the multiphase fluid 101 flows generally upward across the first bend and flows generally downward across the second bend. For example, the multiphase fluid 101 flows generally downward across the first bend and flows generally upward across the second bend. Examples are provided with respect to apparatuses 100E, 100F, 100H, 100I, and 100J shown in FIGS. 1E, 1F, 1H, 1I, and 1J, respectively).

At block 303 f, a third differential pressure (ΔP₃) and a fourth differential pressure (ΔP₄) of the flowing multiphase fluid 101 is measured. The third differential pressure (ΔP₁) is measured in a first direction across a specified differential pressure height (h). The second differential pressure (ΔP₂) is measured in a second direction (different from the first direction) across the specified differential pressure height (h). The first and second directions are such that the multiphase fluid 101 flows generally upward across the specified differential pressure height (h) in one of the first or second directions and flows generally downward across the specified differential pressure height (h) in the other direction. For example, the multiphase fluid 101 flows generally upward across the specified differential pressure height (h) in the first direction and flows generally downward across the specified differential pressure height (h) in the second direction. For example, the multiphase fluid 101 flows generally downward across the specified differential pressure height (h) in the first direction and flows generally upward across the specified differential pressure height (h) in the second direction.

At block 305 f, a frictional pressure loss (ΔP_(f)) of the flowing multiphase fluid 101 is calculated based on the third and fourth differential pressures (ΔP₃, ΔP₄) measured at block 303 f. The frictional pressure loss (ΔP_(f)) can be calculated at block 305 f, for example, by Equation 2h. At block 307 f, an apparent mixture density (ρ_(app)) of the multiphase fluid 101 is estimated based on the third and fourth differential pressures (ΔP₃, ΔP₄) measured at block 303 f. The apparent mixture density (ρ_(app)) can be estimated at block 307 f, for example, by Equation 2d. At block 309 f, a void fraction (λ_(g)) of the multiphase fluid 101 is calculated based on the apparent mixture density (ρ_(app)) estimated at block 307 f. The void fraction (λ_(g)) can be calculated at block 309 f, for example, by Equation 6b.

At block 311 f, various operating conditions and/or properties (such as pressure, temperature, and salinity) about the flowing multiphase fluid 101 are determined. The conditions and/or properties determined at block 311 f can be measured, for example, by external sensors (such as pressure sensors, temperature sensors, and brine samplers). At block 313 f, a water-liquid ratio (R_(WL)) of the flowing multiphase fluid 101 is determined. The water-liquid ratio (R_(WL)) can be determined at block 313 f can be measured, for example, by external sensors (such as dynamic pressure sensors, as in apparatus 100G shown in FIG. 1G). At block 315 f, densities of the various phases (such as density of the gas phase (ρ_(g)), density of the aqueous phase (ρ_(a)), and density of the oil phase (ρ_(o))) of the flowing multiphase fluid 101 are determined. The densities of the various phases of the flowing multiphase fluid 101 can be determined at block 315 f, for example, based on the operating conditions (pressure and temperature) of the flowing multiphase fluid 101 determined at block 311 f. At block 315 f, a density of the liquid phase (ρ_(l), including the aqueous phase and the oil phase) is determined. The density of the liquid phase (ρ_(l)) can be determined at block 315 f based on the densities of the various phases of the flowing multiphase fluid 101 and the water-liquid ratio (R_(WL)) determined at block 313 f. The density of the liquid phase (ρ_(l)) can be determined at block 315 f by Equation 6c. In some implementations, the void fraction (λ_(g)) can be calculated at block 309 f based on the density of the liquid phase (ρ_(l)) and the density of the gas phase (ρ_(g)) determined at block 315 f.

At block 317 f, a gas phase volume fraction (α_(g)) of the multiphase fluid 101 is calculated based on the void fraction (λ_(g)) determined at block 309 f and slip ratio (S) of the flowing multiphase fluid 101. The slip ratio (S) used at block 317 f to determine the gas phase volume fraction (α_(g)) can, for example, be measured by an external sensor or assumed to be substantially equal to 1 as an initial guess. At block 319 f, a mixture density (ρ) of the multiphase fluid 101 is calculated based on the gas phase volume fraction (α_(g)) determined at block 317 f, the density of the gas phase (ρ_(g)) determined at block 315 f, and the density of the liquid phase (ρ_(l)) determined at block 315 f. The mixture density (ρ) can be determined at block 319 f by Equation 6d.

At block 321 f, a friction factor (f_(d)) of the flowing multiphase fluid 101 is determined based on the frictional pressure loss (ΔP_(f)) determined at block 305 f. At block 323 f, a horizontal flow dynamic pressure (q_(h)) and a vertical flow dynamic pressure (q_(v)) are determined. The horizontal flow dynamic pressure (q_(h)) and the vertical flow dynamic pressure (q_(v)) determined at block 323 f can be, for example, measured by external sensors (such as dynamic pressure sensors). In some implementations, the horizontal flow dynamic pressure (q_(h)) is determined by Equation 6f_(h). In some implementations, the vertical flow dynamic pressure (q_(v)) is determined by Equation 6f_(v). In some implementations, the friction factor (f_(d)) is determined at block 321 f by Equation 6g, which depends on the horizontal flow dynamic pressure (q_(h)) and the vertical flow dynamic pressure (q_(v)) determined at block 323 f (for example, by sensors or by Equations 6f_(h) and 6f_(v).

$\begin{matrix} {q_{h} = \frac{{\rho\left( \frac{m_{T}}{\rho A_{h}} \right)}^{2}}{2}} & \left( {6f_{h}} \right) \end{matrix}$ $\begin{matrix} {q_{v} = \frac{{\rho\left( \frac{m_{T}}{\rho A_{v}} \right)}^{2}}{2}} & \left( {6f_{v}} \right) \end{matrix}$ $\begin{matrix} {f_{d} = \frac{\Delta P_{f}}{q_{v}\left( \frac{h}{D_{v}} \right)}} & \left( {6g} \right) \end{matrix}$

In Equation 6f_(h), A_(h) is a cross-sectional flow area for the horizontal flow of the multiphase fluid 101. In Equation 6f_(v), A_(v) is a cross-sectional flow area for the vertical flow of the multiphase fluid 101. In Equation 6g, h is the specified differential pressure height, and D_(v) is the inner diameter of the conduit in which the direction of flow of the multiphase fluid 101 is generally upward or generally downward.

At block 325 f, a total mass flow rate (m_(T)) of the multiphase fluid 101 is determined. The total mass flow rate (m_(T)) can be determined at block 325 f, for example, by a mass flowmeter and/or by any of Equations 4a or 4b. The total mass flow rate (m_(T)) can be determined at block 325 f based on the first and second differential pressures (ΔP₁, ΔP₂) measured at block 301 f, the mixture density (ρ) determined at block 319 f, the friction factor (f_(d)) determined at block 321 f, and the horizontal and vertical flow dynamic pressures (q_(h), q_(v)) determined at block 323 f. In some implementations, the total mass flow rate (m_(T)) is determined at block 325 f by Equation 6h.

$\begin{matrix} {m_{T} = {C_{d}\varepsilon K\sqrt{2{\rho\left( {{{\Delta P_{L}} \pm {\rho{gy}}} - \frac{f_{d}q_{v}y}{D_{v}^{\prime}} - \frac{f_{d}q_{h}x}{D_{h}^{\prime}}} \right)}}}} & \left( {6h} \right) \end{matrix}$

where C_(d) is an empirically characterized discharge coefficient, ε is an expansion factor, K is a fixed geometry factor (K=A/√{square root over (1−β⁴)}, where A is the cross-sectional flow area of the throat of the bend and β is the beta ratio (d/D, where d is the inner diameter downstream of the bend and D is the inner diameter upstream of the bend)), ΔP_(L) is the differential pressure across the bend (selected from one of the first or second differential pressures (ΔP₁ or ΔP₂)), g is the acceleration due to gravity, y is the vertical height difference between the pressure ports for the differential pressure sensor, x is the horizontal distance between the pressure ports for the differential pressure sensor, D′_(h) is the equivalent (effective) inner diameter of the cross-sectional flow area for the horizontal flow of the multiphase fluid 101, and D′_(v) is the equivalent (effective) inner diameter of the cross-sectional flow area for the vertical flow of the multiphase fluid 101.

In some implementations, the blocks in the dotted box of the multiphase fluid flow model 300F (blocks 317 f, 319 f, 321 f, 323 f, 325 f, and determining the slip ratio (S)) are performed simultaneously (for example, in parallel). In some implementations, the blocks in the dotted box of the multiphase fluid flow model 300F (blocks 317 f, 319 f, 321 f, 323 f, 325 f, and determining the slip ratio (S)) are performed iteratively (for example, in series) until the calculated values converge to final values. At block 327 f, various volumetric flow rates of the multiphase fluid 101 are determined. The various volumetric flow rates determined at block 327 f can include a total volumetric flow rate (Q_(T)) of the multiphase fluid 101, a gas phase volumetric flow rate (Q_(g)) of the multiphase fluid 101, a liquid phase volumetric flow rate (Q_(l), including the aqueous phase and the oil phase) of the multiphase fluid 101, an oil phase volumetric flow rate (Q_(o)) of the multiphase fluid 101, and an aqueous phase volumetric flow rate (Q_(a)) of the multiphase fluid 101. The various volumetric flow rates determined at block 327 f can be calculated by Equations 6j_(t), 6j_(g), 6j_(l), 6j_(o), and 6j_(a).

$\begin{matrix} {Q_{T} = \frac{m_{T}}{\rho}} & \left( {6j_{t}} \right) \end{matrix}$ $\begin{matrix} {Q_{g} = {Q_{T}\alpha_{g}}} & \left( {6j_{g}} \right) \end{matrix}$ $\begin{matrix} {Q_{l} = {Q_{T}\left( {1 - \alpha_{g}} \right)}} & \left( {6j_{l}} \right) \end{matrix}$ $\begin{matrix} {Q_{o} = {Q_{T}\left( {1 - R_{WL}} \right)}} & \left( {6j_{o}} \right) \end{matrix}$ $\begin{matrix} {Q_{a} = {Q_{T}R_{WL}}} & \left( {6j_{a}} \right) \end{matrix}$

FIG. 3G is a flowchart depicting an example multiphase fluid flow model 300G. Any one of the apparatuses 100H, 100I, or 100J can, for example, be used to implement the model 300G to characterize a multiphase fluid, such as the multiphase fluid 101. At block 301 g, various operating conditions (such as pressure and temperature) about the flowing multiphase fluid 101 are determined. The various operating conditions determined at block 301 g can be measured, for example, by external sensors, such as pressure sensors and temperature sensors. The various operating conditions determined at block 301 g can include, for example, differential pressure measurements. At block 303 g, a salinity of the aqueous phase of the multiphase fluid 101 is measured. The salinity of the aqueous phase of the multiphase fluid 101 can be measured at block 303 g, for example, using a salinity sensor. In some implementations, the salinity sensor includes electrodes for generating a current to measure conductivity of the aqueous phase of the multiphase fluid 101. At block 305 g, fluid properties (such as density) of the various phases (such as the gas phase, the aqueous phase, and the oil phase) of the multiphase fluid 101 are determined. The fluid properties of the various phases can be determined at block 305 g based on the operating conditions measured at block 301 g and the salinity of the aqueous phase measured at block 303 g. At block 307 g, a void fraction (λ_(g)) of the multiphase fluid 101 is determined. The void fraction (λ_(g)) can be determined at block 307 g, for example, using an external sensor (such as an ultrasonic sensor) and/or by Equation 6b. At block 309 g, a water-liquid ratio (R_(WL)) of the multiphase fluid 101 is determined. The water-liquid ratio (R_(WL)) can be determined at block 309 g, for example, using an external sensor (such as a dynamic pressure sensor). At block 311 g, a density of the liquid phase (ρ_(l), including the aqueous phase and the oil phase) is determined. The density of the liquid phase (ρ_(l)) can be determined at block 311 g based on the densities of the various phases of the flowing multiphase fluid 101 and the water-liquid ratio (R_(WL)) determined at block 309 g. The density of the liquid phase (ρ_(l)) can be determined at block 311 g, for example, by Equation 6c.

At block 313 g, a gas phase volume fraction (α_(g)) of the multiphase fluid 101 is calculated based on the void fraction (λ_(g)) determined at block 307 g and slip ratio (S) of the flowing multiphase fluid 101. The slip ratio (S) used at block 313 g to determine the gas phase volume fraction (α_(g)) can, for example, be measured by an external sensor or assumed to be substantially equal to 1 as an initial guess. At block 315 g, a mixture density (ρ) of the multiphase fluid 101 is calculated based on the gas phase volume fraction (α_(g)) determined at block 313 g, the density of the gas phase (ρ_(g)) determined at block 305 g, and the density of the liquid phase (ρ_(l)) determined at block 311 g. The mixture density (ρ) can be determined at block 315 g, for example, by Equation 6d. At block 317 g, a total mass flow rate (m_(T)) of the multiphase fluid 101 is determined. The total mass flow rate (m_(T)) can be determined at block 317 g, for example, by a mass flowmeter and/or by any of Equations 4a or 4b. In some implementations, the blocks in the dotted box of the multiphase fluid flow model 300G (blocks 313 g, 315 g, 317 g, and determining the slip ratio (S)) are performed simultaneously (for example, in parallel). In some implementations, the blocks in the dotted box of the multiphase fluid flow model 300G (blocks 313 g, 315 g, 317 g, and determining the slip ratio (S)) are performed iteratively (for example, in series) until the calculated values converge to final values.

At block 319 g, various volumetric flow rates of the multiphase fluid 101 are determined. The various volumetric flow rates determined at block 319 g can include a total volumetric flow rate (Q_(T)) of the multiphase fluid 101, a gas phase volumetric flow rate (Q_(g)) of the multiphase fluid 101, a liquid phase volumetric flow rate (Q_(l), including the aqueous phase and the oil phase) of the multiphase fluid 101, an oil phase volumetric flow rate (Q_(o)) of the multiphase fluid 101, and an aqueous phase volumetric flow rate (Q_(a)) of the multiphase fluid 101. The various volumetric flow rates determined at block 319 g can be calculated by Equations 6j_(t), 6j_(g), 6j_(l), 6j_(o), and 6j_(a).

FIG. 4 is a flowchart of an example method 400 for measuring a density of a multiphase fluid 101. Any one of the apparatuses 100A, 100B, 100C, 100D, 100H, 100I, or 100J can be used, for example, to implement the method 400. At block 402, a multiphase fluid (such as the multiphase fluid 101) is flowed through a U-bend (such as the U-bend 102). The multiphase fluid 101 can include a gas phase 101 a, an aqueous phase 101 b′, and an oil phase 101 b″. As described previously, the U-bend 102 includes a first conduit 102 a, a second conduit 102 b, and a connecting conduit 102 c. The first conduit 102 a is parallel to the second conduit 102 b. The connecting conduit 102 c connects the second conduit 102 b to the first conduit 102 a. Each of the first conduit 102 a, the second conduit 102 b, and the connecting conduit 102 c have substantially the same cross-sectional flow area for the multiphase fluid 101 flowing through the U-bend 102. The multiphase fluid 101 flowing through the U-bend 102 flows into the first conduit 102 a, through the connecting conduit 102 c, and out of the second conduit 102 b. At block 404, a first differential pressure (ΔP₁) of the multiphase fluid 101 is measured across a specified differential pressure height (h) along the first conduit 102 a. The first differential pressure (ΔP₁) of the multiphase fluid 101 can be measured at block 404, for example, by the first pair of pressure sensors 104 a, 104 b. At block 406, a second differential pressure (ΔP₂) of the multiphase fluid 101 is measured across the specified differential pressure height (h) along the second conduit 102 b. The second differential pressure (ΔP₂) of the multiphase fluid 101 can be measured at block 406, for example, by the second pair of pressure sensors 104 c, 104 d. At block 408, a mixture density (ρ) of the multiphase fluid 101 is determined at least based on a difference between the first differential pressure (ΔP₁) and the second differential pressure (ΔP₂) of the multiphase fluid 101 (measured at blocks 404 and 406, respectively) and the specified differential pressure height (h). For example, the mixture density (ρ) of the multiphase fluid 101 can be determined at block 406 using Equation 2d.

In some implementations, the method 400 includes determining an aqueous phase volume fraction (α_(a)) of the multiphase fluid 101 at least based on the gas phase volume fraction (α_(g)) of the multiphase fluid 101, the gas density (ρ_(g)) of the gas phase 101 a of the multiphase fluid 101, the aqueous phase density (ρ_(a)) of the aqueous phase 101 b′ of the multiphase fluid 101, the oil phase density (ρ_(o)) of the oil phase 101 b″ of the multiphase fluid 101, and the mixture density (ρ) of the multiphase fluid 101. For example, the gas phase volume fraction (α_(g)), the gas density (ρ_(a)), the aqueous phase density (ρ_(a)), the oil phase density (ρ_(o)), and the mixture density (ρ) of the multiphase fluid 101 can be input to Equation 3b and re-arranged to calculate the aqueous phase volume fraction (α_(a)) of the multiphase fluid 101.

In some implementations, the method 400 includes determining a gas phase volume fraction (α_(g)) of the multiphase fluid 101 at least based on the aqueous phase volume fraction (α_(a)) of the multiphase fluid 101, the gas density (ρ_(g)) of the gas phase 101 a of the multiphase fluid 101, the aqueous phase density (ρ_(a)) of the aqueous phase 101 b′ of the multiphase fluid 101, the oil phase density (ρ_(o)) of the oil phase 101 b″ of the multiphase fluid 101, and the mixture density (ρ) of the multiphase fluid 101. For example, the aqueous phase volume fraction (α_(a)), the gas density (ρ_(g)), the aqueous phase density (ρ_(a)), the oil phase density (ρ_(o)), and the mixture density (ρ) of the multiphase fluid 101 can be input to Equation 3b and re-arranged to calculate the gas phase volume fraction (α_(g)) of the multiphase fluid 101.

FIG. 5 is a flowchart of an example method 500 for measuring a total flow rate of a multiphase fluid 101. Any one of the apparatuses 100E, 100H, 100I, or 100J can be used, for example, to implement the method 500. At block 502, a multiphase fluid (such as the multiphase fluid 101) is flowed from a flow pipe (such as the flow pipe 150) to a U-bend (such as the U-bend 102). The multiphase fluid 101 can include a gas phase 101 a, an aqueous phase 101 b′, and an oil phase 101 b″. As described previously, the U-bend 102 includes a first conduit 102 a, a second conduit 102 b, and a connecting conduit 102 c. The connecting conduit 102 c connects the first conduit 102 a to the second conduit 102 b. The multiphase fluid 101 flowing through the U-bend 102 flows into the first conduit 102 a, through the connecting conduit 102 c, and out of the second conduit 102 b. The flow pipe 150 includes the first pipe 150 a (third conduit) and the second pipe 150 b (fourth conduit). The first pipe 150 a is connected to the first conduit 102 a of the U-bend 102 by a first bend (such as the first angled bend 153 a). The second pipe 150 b is connected to the second conduit 102 b of the U-bend 102 by a second bend (such as the second angled bend 153 b). Flowing the multiphase fluid 101 from the flow pipe 150 to the U-bend 102 at block 502 includes flowing the multiphase fluid 101 from the first pipe 150 a of the flow pipe 150 through the first angled bend 153 a to the first conduit 102 a of the U-bend 102. At block 504, the multiphase fluid 101 is flowed from the U-bend 102 to the flow pipe 150. Flowing the multiphase fluid 101 from the U-bend 102 to the flow pipe 150 at block 504 includes flowing the multiphase fluid 101 from the second conduit 102 b of the U-bend 102 through the second angled bend 153 b to the second pipe 150 b of the flow pipe 150. At block 506, a first differential pressure of the multiphase fluid 101 is measured between a first location on the first pipe 150 a of the flow pipe 150 and a second location on the first conduit 102 a of the U-bend 102. For example, the first differential pressure of the multiphase fluid 101 can be measured at block 506 cooperatively by the first pressure sensor 104 a coupled to the first pipe 150 a of the flow pipe 150 at the first location and by the second pressure sensor 104 b coupled to the first conduit 102 a of the U-bend at the second location. At block 508, a second differential pressure of the multiphase fluid 101 is measured between a third location on the second conduit 102 b of the U-bend 102 and a fourth location on the second pipe 150 b of the flow pipe 150. The second differential pressure of the multiphase fluid 101 can be measured at block 508 cooperatively by the third pressure sensor 104 c coupled to the second conduit 102 b of the U-bend 102 at the third location and by the fourth pressure sensor 104 d coupled to the second pipe 150 b of the flow pipe 150 at the fourth location. At block 510, an overall flow rate of the multiphase fluid 101 is determined at least based on the first differential pressure (measured at block 506), the second differential pressure (measured at block 508), a difference between the first differential pressure and the second differential pressure, and an overall density of the multiphase fluid 101. For example, the total mass flow rate (m_(T1)) of the multiphase fluid 101 can be determined at block 510 using Equation 4a. For example, the total mass flow rate (m_(T2)) of the multiphase fluid 101 can be determined at block 510 using Equation 4b. In some implementations, the method 500 includes comparing the total flow rates (m_(T1), m_(T2)) of the multiphase fluid 101 determined at block 510 (for example, determined by Equations 4a and 4b) to determine an adjustment of the discharge coefficients (for example, C_(d1), C_(d2)).

FIG. 6 is a flowchart of an example method 600 for measuring a mixture density, a total flow rate, and phase flow rates of a multiphase fluid 101. Any of the apparatuses 100H, 100I, or 100J can be used, for example, to implement the method 600. At block 602, a multiphase fluid (such as the multiphase fluid 101) is flowed from a flow pipe (such as the flow pipe 150) to a U-bend (such as the U-bend 102). The multiphase fluid 101 can include a gas phase 101 a, an aqueous phase 101 b′, and an oil phase 101 b″. As described previously, the U-bend 102 includes a first conduit 102 a, a second conduit 102 b, and a connecting conduit 102 c. The connecting conduit 102 c connects the first conduit 102 a to the second conduit 102 b. The multiphase fluid 101 flowing through the U-bend 102 flows into the first conduit 102 a, through the connecting conduit 102 c, and out of the second conduit 102 b. The flow pipe 150 includes the first pipe 150 a (third conduit) and the second pipe 150 b (fourth conduit). The first pipe 150 a is connected to the first conduit 102 a of the U-bend 102 by a first bend (such as the first angled bend 153 a). The second pipe 150 b is connected to the second conduit 102 b of the U-bend 102 by a second bend (such as the second angled bend 153 b). Flowing the multiphase fluid 101 from the flow pipe 150 to the U-bend 102 at block 602 includes flowing the multiphase fluid 101 from the first pipe 150 a of the flow pipe 150 through the first angled bend 153 a to the first conduit 102 a of the U-bend 102. At block 604, the multiphase fluid 101 is flowed from the U-bend 102 to the flow pipe 150. Flowing the multiphase fluid 101 from the U-bend 102 to the flow pipe 150 at block 604 includes flowing the multiphase fluid 101 from the second conduit 102 b of the U-bend 102 through the second angled bend 153 b to the second pipe 150 b of the flow pipe 150. At block 606, a first differential pressure of the multiphase fluid 101 is measured between a first location on the first pipe 150 a of the flow pipe 150 and a second location on the first conduit 102 a of the U-bend 102. For example, the first differential pressure of the multiphase fluid 101 can be measured at block 606 cooperatively by the first pressure sensor 104 a coupled to the first pipe 150 a of the flow pipe 150 at the first location and by the second pressure sensor 104 b coupled to the first conduit 102 a of the U-bend at the second location. At block 608, a second differential pressure of the multiphase fluid 101 is measured between a third location on the second conduit 102 b of the U-bend 102 and a fourth location on the second pipe 150 b of the flow pipe 150. The second differential pressure of the multiphase fluid 101 can be measured at block 608 cooperatively by the third pressure sensor 104 c coupled to the second conduit 102 b of the U-bend 102 at the third location and by the fourth pressure sensor 104 d coupled to the second pipe 150 b of the flow pipe 150 at the fourth location. At block 610, a third differential pressure of the multiphase fluid 101 is measured across a specified height (such as the specified height, h) along the longitudinal length of the first conduit 102 a of the U-bend 102. The third differential pressure of the multiphase fluid 101 can be measured at block 610 by the third pair of pressure sensors 104 e, 104 f coupled to the first conduit 102 a. At block 612, a fourth differential pressure of the multiphase fluid 101 is measured across the specified height (h) along the longitudinal length of the second conduit 102 b of the U-bend 102. The fourth differential pressure of the multiphase fluid 101 can be measured at block 612 by the fourth pair of pressure sensors 104 g, 104 h coupled to the second conduit 102 b. At block 614, a mixture density (ρ) of the multiphase fluid 101 is determined at least based on a difference between the third differential pressure (ΔP₃) and the fourth differential pressure (ΔP₄) of the multiphase fluid 101 (measured at blocks 610 and 612, respectively) and the specified height (h). For example, the mixture density (ρ) of the multiphase fluid 101 can be determined at block 614 using Equation 2d, where ΔP₁ is substituted for the third differential pressure (ΔP₃) measured by the third pair of pressure sensors 104 e, 104 f at block 610, ΔP₂ is substituted for the fourth differential pressure (ΔP₄) measured by the fourth pair of pressure sensors 104 g, 104 h at block 612, and h₂ is substituted for the specified height (h). At block 616, a total flow rate of the multiphase fluid 101 is determined at least based on the first differential pressure (measured at block 606), the second differential pressure (measured at block 608), and a mixture density of the multiphase fluid 101 (determined at block 614). For example, the total mass flow rate (m_(T1)) of the multiphase fluid 101 can be determined at block 616 using Equation 4a. For example, the total mass flow rate (m_(T2)) of the multiphase fluid 101 can be determined at block 616 using Equation 4b.

In some implementations, the method 600 includes comparing the total flow rates (m_(T1), m_(T2)) of the multiphase fluid 101 determined at block 616 (for example, determined by Equations 4a and 4b) to determine an adjustment of the discharge coefficients (C_(d1), C_(d2)). In implementation in which the method 600 is performed by the apparatus 100G, the method 600 can include measuring various dynamic pressures of the multiphase fluid 101 flowing through the U-bend 102. The dynamic pressures of the multiphase fluid 101 can be measured, for example, by the dynamic pressure sensors 106 a, 106 b, 106 c. The dynamic pressure measurements taken by the dynamic pressure sensors 106 a, 106 b, 106 c can be cross-correlated (for example, by the processor 705) to determine a speed of sound of the multiphase fluid 101 flowing through the U-bend 102. The speed of sound of the multiphase fluid 101 and the mixture density (ρ) of the multiphase fluid 101 can then be used to determine the water-liquid ratio of the multiphase fluid 101. For example, the speed of sound of the multiphase fluid 101 and the mixture density (ρ) of the multiphase fluid 101 can be compared to a plotted graph of various phase curves (gas/oil, gas/aqueous, oil/aqueous) to determine the water-liquid ratio.

FIG. 7 is a block diagram of an example computer 700 used to provide computational functionalities associated with described algorithms, methods, functions, processes, flows, and procedures, as described in this specification, according to an implementation. The illustrated computer 700 is intended to encompass any computing device such as a server, desktop computer, laptop/notebook computer, one or more processors within these devices, or any other processing device, including physical or virtual instances (or both) of the computing device. Additionally, the computer 700 can include a computer that includes an input device, such as a keypad, keyboard, touch screen, or other device that can accept user information, and an output device that conveys information associated with the operation of the computer 700, including digital data, visual, audio information, or a combination of information.

The computer 700 includes an interface 704. Although illustrated as a single interface 704 in FIG. 7 , two or more interfaces 704 may be used according to particular needs, desires, or particular implementations of the computer 700. Although not shown in FIG. 7 , the computer 700 can be communicably coupled with a network. The interface 704 is used by the computer 700 for communicating with other systems that are connected to the network in a distributed environment. Generally, the interface 704 comprises logic encoded in software or hardware (or a combination of software and hardware) and is operable to communicate with the network. More specifically, the interface 704 may comprise software supporting one or more communication protocols associated with communications such that the network or interface's hardware is operable to communicate physical signals within and outside of the illustrated computer 700. The interface 704 can include a control interface, which can be used to couple the computer 700 to controls. In some implementations, the control interface is a bank of relays, a bank of MOSFET power controllers, a serial peripheral interface (SPI), or a Fieldbus, and the like. The interface 704 can include a sensor interface, which can be used to couple the computer 700 to sensors. In some implementations, the sensor interface is a bank of analog-to-digital converters (ADCs), and I2C bus, a serial peripheral interface (SPI) bus, or a Fieldbus, and the like. The interface 704 can include a human machine interface, which can be used by a user to interact with the computer 700. In some implementations, the human machine interface includes a monitor or a touch screen that is configured to display information, for example, to a user.

The computer 700 includes a processor 705. The processor 705 may be a microprocessor, a multi-core processor, a multithreaded processor, an ultra-low-voltage processor, an embedded processor, or a virtual processor. In some embodiments, the processor 705 may be part of a system-on-a-chip (SoC) in which the processor 705 and the other components of the computer 700 are formed into a single integrated electronics package. In some implementations, the processor 705 may include processors from Intel® Corporation of Santa Clara, California, from Advanced Micro Devices, Inc. (AMD) of Sunnyvale, California, or from ARM Holdings, LTD., Of Cambridge, England. Any number of other processors from other suppliers may also be used. Although illustrated as a single processor 705 in FIG. 7 , two or more processors may be used according to particular needs, desires, or particular implementations of the computer 700. Generally, the processor 705 executes instructions and manipulates data to perform the operations of the computer 700 and any algorithms, methods, functions, processes, flows, and procedures as described in this specification. The processor 705 may communicate with other components of the computer 700 over a bus. The bus may include any number of technologies, such as industry standard architecture (ISA), extended ISA (EISA), peripheral component interconnect (PCI), peripheral component interconnect extended (PCIx), PCI express (PCIe), or any number of other technologies. The bus may be a proprietary bus, for example, used in an SoC based system. Other bus technologies may be used, in addition to, or instead of, the technologies above.

The computer 700 can also include a database 706 that can hold data for the computer 700 or other components (or a combination of both) that can be connected to the network. Although illustrated as a single database 706 in FIG. 7 , two or more databases (of the same or combination of types) can be used according to particular needs, desires, or particular implementations of the computer 700 and the described functionality. While database 706 is illustrated as an integral component of the computer 700, database 706 can be external to the computer 700. The database 706 can be used for the persistent storage of information, such as data, applications, operating systems, and so forth. The database 706 may be a nonvolatile RAM, a solid-state disk drive, or a flash drive, among others. In some implementations, the database 306 will include a hard disk drive, such as a micro hard disk drive, a regular hard disk drive, or an array of hard disk drives, for example, associated with a DCS or a cloud server.

The computer 700 also includes a memory 707 that can hold data for the computer 700 or other components (or a combination of both) that can be connected to the network. Although illustrated as a single memory 707 in FIG. 7 , two or more memories 707 (of the same or combination of types) can be used according to particular needs, desires, or particular implementations of the computer 700 and the described functionality. While memory 707 is illustrated as an integral component of the computer 700, memory 707 can be external to the computer 700. The memory 707 can be a transitory or non-transitory storage medium. In some implementations, such as in PLCs and other process control units, the memory 707 is integrated with the database 706 used for long-term storage of programs and data. The memory 707 can include any number of volatile and nonvolatile memory devices, such as volatile random-access memory (RAM), static random-access memory (SRAM), flash memory, and the like. In smaller devices, such as PLCs, the memory 707 may include registers associated with the processor 705 itself.

The memory 707 stores computer-readable instructions executable by the processor 705 that, when executed, cause the processor 705 to perform operations, such as receiving a first differential pressure signal from the first pair of pressure sensors 104 a, 104 b; receiving a second differential pressure signal from the second pair of pressure sensors 104 c, 104 d; and determining an overall density of the multiphase fluid 101 at least based on a difference between the first differential pressure and the second differential pressure of the multiphase fluid 101 and the specified differential pressure height, h. The computer 700 can also include a power supply 714. The power supply 714 can include a rechargeable or non-rechargeable battery that can be configured to be either user- or non-user-replaceable. The power supply 714 can be hard-wired. There may be any number of computers 700 associated with, or external to, a computer system containing computer 700, each computer 700 communicating over the network. Further, the term “client,” “user,” “operator,” and other appropriate terminology may be used interchangeably, as appropriate, without departing from this specification. Moreover, this specification contemplates that many users may use one computer 700, or that one user may use multiple computers 700.

While this specification contains many specific implementation details, these should not be construed as limitations on the scope of what may be claimed, but rather as descriptions of features that may be specific to particular implementations. Certain features that are described in this specification in the context of separate implementations can also be implemented, in combination, in a single implementation. Conversely, various features that are described in the context of a single implementation can also be implemented in multiple implementations, separately, or in any sub-combination. Moreover, although previously described features may be described as acting in certain combinations and even initially claimed as such, one or more features from a claimed combination can, in some cases, be excised from the combination, and the claimed combination may be directed to a sub-combination or variation of a sub-combination.

As used in this disclosure, the terms “a,” “an,” or “the” are used to include one or more than one unless the context clearly dictates otherwise. The term “or” is used to refer to a nonexclusive “or” unless otherwise indicated. The statement “at least one of A and B” has the same meaning as “A, B, or A and B.” In addition, it is to be understood that the phraseology or terminology employed in this disclosure, and not otherwise defined, is for the purpose of description only and not of limitation. Any use of section headings is intended to aid reading of the document and is not to be interpreted as limiting; information that is relevant to a section heading may occur within or outside of that particular section.

As used in this disclosure, the term “about” or “approximately” can allow for a degree of variability in a value or range, for example, within 10%, within 5%, or within 1% of a stated value or of a stated limit of a range.

As used in this disclosure, the term “substantially” refers to a majority of, or mostly, as in at least about 50%, 60%, 70%, 80%, 90%, 95%, 96%, 97%, 98%, 99%, 99.5%, 99.9%, 99.99%, or at least about 99.999% or more.

Values expressed in a range format should be interpreted in a flexible manner to include not only the numerical values explicitly recited as the limits of the range, but also to include all the individual numerical values or sub-ranges encompassed within that range as if each numerical value and sub-range is explicitly recited. For example, a range of “0.1% to about 5%” or “0.1% to 5%” should be interpreted to include about 0.1% to about 5%, as well as the individual values (for example, 1%, 2%, 3%, and 4%) and the sub-ranges (for example, 0.1% to 0.5%, 1.1% to 2.2%, 3.3% to 4.4%) within the indicated range. The statement “X to Y” has the same meaning as “about X to about Y,” unless indicated otherwise. Likewise, the statement “X, Y, or Z” has the same meaning as “about X, about Y, or about Z,” unless indicated otherwise.

Particular implementations of the subject matter have been described. Other implementations, alterations, and permutations of the described implementations are within the scope of the following claims as will be apparent to those skilled in the art. While operations are depicted in the drawings or claims in a particular order, this should not be understood as requiring that such operations be performed in the particular order shown or in sequential order, or that all illustrated operations be performed (some operations may be considered optional), to achieve desirable results. In certain circumstances, multitasking or parallel processing (or a combination of multitasking and parallel processing) may be advantageous and performed as deemed appropriate.

Moreover, the separation or integration of various system modules and components in the previously described implementations should not be understood as requiring such separation or integration in all implementations, and it should be understood that the described components and systems can generally be integrated together or packaged into multiple products.

Accordingly, the previously described example implementations do not define or constrain the present disclosure. Other changes, substitutions, and alterations are also possible without departing from the spirit and scope of the present disclosure. 

What is claimed is:
 1. A system comprising: a U-bend configured to flow a multiphase fluid, the U-bend comprising: a first conduit; a second conduit, wherein a cross-sectional flow area of the first conduit is substantially the same as a cross-sectional flow area of the second conduit; and a connecting conduit connecting the first conduit to the second conduit; a third conduit configured to flow the multiphase fluid, the third conduit having a cross-sectional flow area that is different from the cross-sectional flow area of the first conduit; a first bend connecting the third conduit to the first conduit; a fourth conduit configured to flow the multiphase fluid, the fourth conduit parallel to and in-line with the third conduit, the fourth conduit having a cross-sectional flow area that is substantially the same as the cross-sectional flow area of the third conduit; a second bend connecting the fourth conduit to the second conduit; a first differential pressure sensor coupled to the third conduit and to the first conduit, wherein the first differential pressure sensor is configured to measure a first differential pressure of the multiphase fluid across the first bend; a second differential pressure sensor coupled to the second conduit and to the fourth conduit, wherein the second differential pressure sensor is configured to measure a second differential pressure of the multiphase fluid across the second bend; and a computer, comprising: a processor communicatively coupled to the first differential pressure sensor and to the second differential pressure sensor; and a computer-readable storage medium coupled to the processor and storing programming instructions for execution by the processor, the programming instructions instructing the processor to perform operations comprising: receiving a first differential pressure signal from the first differential pressure sensor, the first differential pressure signal representing the first differential pressure of the multiphase fluid; receiving a second differential pressure signal from the second differential pressure sensor, the second differential pressure signal representing the second differential pressure of the multiphase fluid; receiving, as input, a mixture density of the multiphase fluid; and determining a total flow rate of the multiphase fluid at least based on the first differential pressure, the second differential pressure, and the mixture density of the multiphase fluid.
 2. The system of claim 1, wherein the first differential pressure sensor is coupled to the third conduit at a first location and to the first conduit at a second location, and the first location and the second location are separated by a specified height with respect to gravity.
 3. The system of claim 2, wherein the second differential pressure sensor is coupled to the second conduit at a third location and to the fourth conduit at a fourth location, and the third location and the second location are separated by the specified height with respect to gravity.
 4. The system of claim 3, wherein the total flow rate of the multiphase fluid is determined as a first total mass flow rate of the multiphase fluid at least based on a first differential pressure equation: m _(T1) =C _(d1)×ε₁ ×K ₁×√{square root over (2×ρ×(ΔP ₁ −ρ×g×h))}, wherein m_(T1) is the first total mass flow rate of the multiphase fluid, C_(d1) is a first discharge coefficient, ε₁ is a first expansion factor, K₁ is a first fixed geometry factor proportional to the cross-sectional flow area of the first conduit, ΔP₁ is the first differential pressure, ρ is the mixture density of the multiphase fluid, g is an acceleration due to gravity, and h is the specified height.
 5. The system of claim 4, wherein the total flow rate of the multiphase fluid is determined as a second total mass flow rate of the multiphase fluid at least based on a second differential pressure equation: m _(T2) =C _(d2)×ε₂ ×K ₂×√{square root over (2×ρ×(ΔP ₂ −ρ×g×h))}, wherein m_(T2) is the second total mass flow rate of the multiphase fluid, C_(d2) is a second discharge coefficient, ε₂ is a second expansion factor, K₂ is a second fixed geometry factor proportional to the cross-sectional flow area of the fourth conduit, ΔP₂ is the second differential pressure, ρ is the mixture density of the multiphase fluid, g is the acceleration due to gravity, and h is the specified height.
 6. The system of claim 5, wherein the operations performed by the processor comprise comparing the first total mass flow rate of the multiphase fluid and the second total mass flow rate of the multiphase fluid to determine an adjustment of the first discharge coefficient (C_(d1)) or an adjustment of the second discharge coefficient (C_(d2)).
 7. The system of claim 6, wherein the first total mass flow rate (m_(T1)) is equal to the second total mass flow rate (m_(T2)), the mixture density of the multiphase fluid is a first mixture density of the multiphase fluid, and the operations performed by the processor comprise determining a second mixture density of the multiphase fluid at least based on a combined differential pressure equation: ${\rho_{2} = \frac{{Z \times \Delta P_{1}} - {\Delta P_{2}}}{g \times \left( {1 + Z} \right)}},$ wherein ρ₂ is the second mixture density of the multiphase fluid, ΔP₁ is the first differential pressure, ΔP₂ is the second differential pressure, g is the acceleration due to gravity, and Z is defined by: ${Z = \left( \frac{C_{d1} \times \varepsilon_{1} \times K_{1}}{C_{d2} \times \varepsilon_{2} \times K_{2}} \right)^{2}},$ wherein C_(d1) is the first discharge coefficient, ε₁ is the first expansion factor, K₁ is the first fixed geometry factor, C_(d2) is the second discharge coefficient, ε₂ is the second expansion factor, and K₂ is the second fixed geometry factor.
 8. The system of claim 7, wherein the first differential pressure sensor is a first circumferential pressure sensor that spans an entire circumference of the first conduit.
 9. The system of claim 8, wherein the second differential pressure sensor is a second circumferential pressure sensor that spans an entire circumference of the second conduit.
 10. The system of claim 9, wherein each of the first bend and the second bend are 90-degree angled bends.
 11. A method comprising: flowing a multiphase fluid from a flow pipe to a U-bend, wherein the multiphase fluid comprises a gas phase, an aqueous phase, and an oil phase, and the U-bend comprises: a first conduit; a second conduit; and a connecting conduit connecting the first conduit to the second conduit, wherein the multiphase fluid flowing through the U-bend flows into the first conduit, through the connecting conduit, and out of the second conduit; flowing the multiphase fluid from the U-bend to the flow pipe, wherein: the flow pipe comprises a third conduit and a fourth conduit; the third conduit of the flow pipe is connected to the first conduit of the U-bend by a first bend; the fourth conduit of the flow pipe is connected to the second conduit of the U-bend by a second bend; flowing the multiphase fluid from the flow pipe to the U-bend comprises flowing the multiphase fluid from the third conduit of the flow pipe through the first bend to the first conduit of the U-bend; and flowing the multiphase fluid from the U-bend to the flow pipe comprises flowing the multiphase fluid from the second conduit of the U-bend through the second bend to the fourth conduit of the flow pipe; measuring a first differential pressure of the multiphase fluid between a first location on the third conduit of the flow pipe and a second location on the first conduit of the U-bend, wherein the first location and the second location are separated by a specified height along the first conduit of the U-bend; measuring a second differential pressure of the multiphase fluid between a third location on the second conduit of the U-bend and a fourth location on the fourth conduit of the flow pipe, wherein the third location and the fourth location are separated by the specified height along the second conduit of the U-bend; and determining a total flow rate of the multiphase fluid at least based on the first differential pressure, the second differential pressure, and a mixture density of the multiphase fluid.
 12. The method of claim 11, wherein the total flow rate of the multiphase fluid is determined as a first total mass flow rate of the multiphase fluid at least based on a first differential pressure equation: m _(T1) =C _(d1)×ε₁ ×K ₁×√{square root over (2×ρ×(ΔP ₁ −ρ×g×h))}, wherein m_(T1) is the first total mass flow rate of the multiphase fluid, C_(d1) is a first discharge coefficient, ε₁ is a first expansion factor, K₁ is a first fixed geometry factor proportional to the cross-sectional flow area of the first conduit, ΔP₁ is the first differential pressure, ρ is the mixture density of the multiphase fluid, g is an acceleration due to gravity, and h is the specified height.
 13. The method of claim 12, wherein the total flow rate of the multiphase fluid is determined as a second total mass flow rate of the multiphase fluid at least based on a second differential pressure equation: m _(T2) =C _(d2)×ε₂ ×K ₂×√{square root over (2×ρ×(ΔP ₂ −ρ×g×h))}, wherein m_(T2) is the second total mass flow rate of the multiphase fluid, C_(d2) is a second discharge coefficient, ε₂ is a second expansion factor, K₂ is a second fixed geometry factor proportional to the cross-sectional flow area of the fourth conduit, ΔP₂ is the second differential pressure, ρ is the mixture density of the multiphase fluid, g is the acceleration due to gravity, and h is the specified height.
 14. The method of claim 13, wherein the first total mass flow rate (m_(T1)) is equal to the second total mass flow rate (m_(T2)), the mixture density of the multiphase fluid is a first mixture density of the multiphase fluid, and the operations performed by the processor comprise determining a second mixture density of the multiphase fluid at least based on a combined differential pressure equation: ${\rho_{2} = \frac{{Z \times \Delta P_{1}} - {\Delta P_{2}}}{g \times \left( {1 + Z} \right)}},$ wherein ρ₂ is the second mixture density of the multiphase fluid, ΔP₁ is the first differential pressure, ΔP₂ is the second differential pressure, g is the acceleration due to gravity, and Z is defined by: ${Z = \left( \frac{C_{d1} \times \varepsilon_{1} \times K_{1}}{C_{d2} \times \varepsilon_{2} \times K_{2}} \right)^{2}},$ wherein C_(d1) is the first discharge coefficient, ε₁ is the first expansion factor, K₁ is the first fixed geometry factor, C_(d2) is the second discharge coefficient, ε₂ is the second expansion factor, and K₂ is the second fixed geometry factor.
 15. The method of claim 14, comprising comparing the first total mass flow rate of the multiphase fluid and the second total mass flow rate of the multiphase fluid to determine an adjustment of the first discharge coefficient (C_(d1)) or an adjustment of the second discharge coefficient (C_(d2)).
 16. The method of claim 15, wherein the first differential pressure is measured by a first differential pressure sensor coupled to the third conduit of the flow pipe at the first location and coupled to the first conduit of the U-bend at the second location, the second differential pressure is measured by a second differential pressure sensor coupled to the second conduit of the U-bend at the third location and coupled to the fourth conduit of the flow pipe at the fourth location.
 17. The method of claim 16, wherein the first differential pressure sensor is a first circumferential pressure sensor that spans an entire circumference of the first conduit.
 18. The method of claim 17, wherein the second differential pressure sensor is a second circumferential pressure sensor that spans an entire circumference of the second conduit.
 19. The method of claim 18, wherein each of the first bend and the second bend are 90-degree angled bends. 